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Chapter 1 



Are you entangled? 



Introduction 

1.1 Overview and Objectives 

In the last two decades quantum information theory has emerged as a branch 
of quantum physics which links concepts of quantum physics to ideas from in- 
formation theory. On the one hand it became evident that the language from 
classical information theory is appropriate to grasp several important open ques- 
tions in quantum mechanics. On the other hand, it was recognised that quan- 
tum mechanical systems would allow for processing of classical information that 
could not be realised with physical systems belonging to the realm of classical 
physics. Most prominently perhaps, factoring of large numbers is a problem 
that can be solved in polynomial time on a quantum computer. In computer 
science, this problem belongs to the complexity class of NP problems, which 
means that there is no known polynomial time algorithm for a classical com- 
puter. Frankly, to factor large numbers is simply not possible in any reasonable 
time scale with a classical computer. The insight that there exist quantum 
algorithms that so massively outperform the best known classical algorithms 
resulted in a strong interest in quantum computation. Soon other subjects of 
information theory were reviewed from the quantum physics perspective, and 
many other applications were found. 

Furthermore it was recognised that entanglement plays a major role for those 
applications which can not be realised with classical schemes or setups. Today, 
entanglement is seen as a resource, which is consumed when using it for quantum 
information processing. The dictionary describes entangled as follows: 

If something is entangled in something else, such as a rope, wire, or 
net, it is caught in it very firmly. 

If you are entangled in something, you are involved in difficulties 
from which it is hard to escape. 

If you are entangled with someone, you are involved in a relationship 
with them that causes problems or difficulties. 

In quantum physics, entangled states can be understood as states where degrees 
of freedom of parts of the system are tightly correlated. In contrast to clas- 
sical correlations, entanglement can not be generated by a preparation using 
only local operations on parts of the system and thus offers new dimensions 
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for information theory and leads to manifestly counterintuitive consequences. 
Since entanglement is the necessary ingredient for so many applications, one 
naturally wants to know how to quantify it. Different entanglement measures 
were proposed, and the concepts of distillable entanglement and entanglement 
cost were introduced. 

Until recently, research in the field of the theory of entanglement and its appli- 
cations has almost exclusively dealt with finite dimensional quantum systems 
such as two- level systems, in analogy to the classical bits called qubits. The 
theoretical concept of two-level systems can be realised in nature by the two 
directions of a spin in a magnetic field or by the polarisation degrees of freedom 
of light. It was, for example, shown that teleportation of an unknown two-level 
state between two parties is possible, if both share a pair of entangled states. 
The experimental realisation was done in 1998, Ref. [1] . 

Discrete systems, like the two-level systems, are associated with finite-dimens- 
ional Hilbert spaces. In recent years, continuous variable systems (CVS) came 
into the focus of research and fruitful attempts were done to employ the con- 
cepts and insights from discrete systems on continuous variable states, living 
in an infinite dimensional Hilbert space. Examples are the quantised electro- 
magnetic field, known as a set of harmonic oscillators, the vibrational degree 
of freedom of ions in a trap, or mechanical oscillators in the quantum domain. 
The simplest nontrivial yet very important class of CVS are the Gaussian states 
determined by their first and second moments alone. They show most of the 
interesting effects like entanglement and squeezing and are mathematically easy 
to handle. Furthermore Gaussian states are experimentally available in optical 
settings. 

We will investigate the entanglement properties of Gaussian states and formu- 
late them on covariance matrices, associated to the states. Furthermore, we 
do not only want to quantify the degree of entanglement, but also to estimate 
how much a given state is entangled. This question is not easy to answer even 
in the case when the state is fully known, but additionally one is interested 
in good strategies of how to estimate or give a minimum amount of entangle- 
ment of an unknown or partly known state. The young theory of entanglement 
witnesses is introduced and provides a powerful tool to answer this question. 
As we will present in the sixth chapter, it is indeed possible to use entangle- 
ment witnesses in order to estimate the degree of entanglement of an unknown 
state with experimentally available measurements. To quantify entanglement 
we will throughout this thesis use the logarithmic negativity, being so far the 
only computable entanglement measure for general states. 



1.2 Outline 



This thesis is divided in seven chapters, and an appendix collecting some small 
definitions and proofs. The material is organised as follows: 
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Basics 

In this chapter several concepts of continuous variable quantum mechanics are 
reviewed, namely the commutation relations, the concept of phase space and of 
canonical (symplectic) transformations. We introduce the moments of a state 
and the quantum analogue to the classical characteristic function. Finally, we 
define entanglement and state the famous PPT-criterion. 

Gaussian States 

In the third chapter we introduce the set of Gaussian states starting with some 
familiar examples as the coherent states of a harmonic oscillator. We define the 
covariance matrix as the major quantity determining most of the important 
properties of a Gaussian state. Furthermore two useful normal forms of covari- 
ance matrices will be introduced connected to invariant quantities, the sym- 
plectic eigenvalues and the Simon invariants, whose properties we will analyse. 
Additionally we see that the Gaussian states maximise the entropy. 

Gaussian Operations 

This chapter mainly focuses on homodyne measurements but also shortly dis- 
cusses the allowed operations one could perform on a Gaussian state without 
destroying this property. We show that there exists an efficient setup to project 
an arbitrary state on a coherent state using homodyne measurements and other 
Gaussian operations. As we see in the subsequent chapter, the derived setup is 
necessary for our scheme to determine the degree of entanglement. 

Measuring Entanglement of Two-Mode States 

In this and the subsequent chapter the methods of estimating the degree of 
entanglement are introduced. First, we present a scheme of how to measure the 
degree of entanglement of an unknown quantum state. We propose a setup to 
give the exact logarithmic negativity of a Gaussian two-mode state with only 
nine kinds of measurements instead of estimating all entries of the covariance 
matrix with homodyne detections. 

Entanglement Witnesses 

In this chapter we introduce the theory of entanglement witnesses and show 
that it is a helpful tool to answer our question for entanglement properties 
of states. We notice that the required measurements are easily available in 
the laboratory and show that it is possible not only to decide experimentally 
whether an unknown given state is entangled or not, but also to give a lower 
bound of its logarithmic negativity. To do this we introduce the new quanti- 
ties p— separability and minimal entanglement witness, which characterise the 
entanglement properties of covariance matrices. Finally, we present a general 
strategy to estimate the degree of entanglement of an unknown Gaussian state, 
using the concept of entanglement witnesses. As the closing words, we give an 
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example and shortly discuss the Duan-criterion as a special case of the general 
criterion we stated. 

Summary and Outlook 

We finally summarise the main results of this thesis and give some possible 
issues of further research. 



Chapter 2 



Basics 



In this thesis we will discuss aspects of quantum information theory (QIT) 
with continuous variable states - especially Gaussian states. These states can 
be described in a different formalism than the one of the familiar discrete qubit 
systems. We therefore start off with some basic concepts for continuous variable 
states (CVS). 



2.1 Commutation Relations 

In classical mechanics a system is described completely by its degrees of freedom, 
usually the coordinates (x n ) and momenta (p n ) of all N particles (with n = 
1, N in one-dimensional space). These canonical variables fulfil the canonical 
relations 



{Ri, Rj} — Oi 



(2.1) 



where { , } denotes the Poisson bracket defined in Def. A.l, with R E R 2Ar 
being the vector of canonical variables R = (x±,p±, ... ,xn,pn) t and a the 
symplectic matrix defined as: 



a 



N 

e 

i=i 



o 1 
-1 o 



1 
-1 



V 








\ 




1 
-1 0/ 



(2.2) 



The phase space (H 2N , a) spanned by the canonical coordinates is then a sym- 
plectic vector space (see Def. A. 2) with a symplectic scalar product 



which is preserved under symplectic basis transformations, defined in the next 
section. All a appearing in this thesis denote symplectic matrices of appropri- 
ate size. 
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In quantum mechanics we describe the canonical degrees of freedom by the set 
of operators {x n ,p n } n or equivalently by the set of ladder operators {a n ,a n } n 
with n = 1, . . . , N. The generalised position operators x n act in the position 
representation on a vector tp G TC as multiplications 

x n tp(x) = x n tp(x) 

while the generalised momentum operators p n act as derivatives 

. dip(x) 



Pn tp(x) 



dx n 



These are unbounded linear, selfadjoint operators on the Hilbert space TL = 
£ 2 (H N , C). When quantising a system the Poisson brackets translate to com- 
mutators and the classical canonical relations become the canonical commuta- 
tion relations (CCRs): 

[R j ,R k ]=ia jk l (2.3) 

with R = (R\, ... ,R,2n) T = (^i)Plj ••• :Xn,Pn) T and a being the symplectic 
matrix. Sometimes we will also use the ladder operators, connected to the 
position and momentum operators via 

a n := — -j=— and < := — -j=— , (2.4) 
x n = and p n = , (2.5) 

and the translated CCRs read a k ] = 5jk 1, with all other commutators equal 
to zero. 



2.2 Symplectic Transformations 

When changing our coordinate system or doing any transformation on a quan- 
tum mechanical state we use canonical transformations, T : R R' , which 
leave the basic kinematic relations unchanged, e.g., they leave the CCRs un- 
changed: [R'j,R' k ] = icrjkt. Apart from shifts of the coordinate system's origin 
the simplest canonical transformations are linear homogeneous ones. These are 
called symplectic transformations (see Ref. [2, 3, 4]) which form the real sym- 
plectic group, being one of the three major semisimple Lie groups besides the 
real orthogonal groups and the complex unitary groups. 

Definition 2.1 The symplectic group will be denoted by 

Sp(2N, R) = {S\ real 2N x 2N matrix , SaS T = a}, (2.6) 

where the last condition comes from the CCRs, which shall still be fulfilled 
when expressing the new coordinates as linear functions of the old ones via 



S : R^ R' = SR, 



(2.7) 
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with S being a symplectic transformation. 

We find the following properties for the matrices S G Sp(2N, R) : 

1. dim(Sp(2N, R)) = N(2N + 1), 

2. a G Sp(2N, R), (2.8) 

3. 5 G 5p(2iV, R) => -5, 5 -1 , 5 T G 5p(2JV, R), 

4. detS = +l. 

The first property comes from the restrictions on the entries of S given by 
SaS T = a, 2 and 3 are easy to check while the last condition is rather subtle 
and will be proved in the appendix (see Def. A. 2 and the following discussion). 
Note, that all 2 x 2 matrices with determinant one are symplectic matrices, 
especially all matrices from SO(2, R). 



2.2.1 Subsets of the Symplectic Group 

Sometimes it is useful to divide a symplectic transformation in transformations 
on subsystems. With S = ( ^ ^ J arm c in a basis where a looks like 



IjVxJV 

■Ijvxjv 



, the properties of the symplectic S may be translated 



to: 



AB T , CD T , A T C, B T D symmetric, 

AD T - BC T = l NxN , (2.9) 
A T D - C T B = t NxN . 

Obviously, the symplectic transformations do in general not preserve angles 
between vectors. We therefore introduce the important compact, maximally 
complete subgroup K{N) = Sp{2N, R) n SO(2N) of Sp(2N, R) and find that 
K(N) is isomorphic to the unitary group U(N). 

That is because every unitary NxN matrix can be expanded U = X + iY 
where X and Y are real NxN matrices with 

X T X + Y T Y =t NxN = XX T + YY T , 
X T Y, XY T symmetric. 

The submatrices of any S G K{N) fulfil exactly the above conditions required 
for any element of U(N). Hence K(N) is isomorphic to U(N). We write the 
set K(N) as : 

k(n) = \k\k =( y x ) with x + iYe c/ ( iV ); X ' y real }• ( 2 - 10 ) 



These linear homogeneous transformations will be called passive transforma- 
tions because they conserve the energy of the system. 
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Example 1 : The two-mode (50:50) beam splitter is an element 
of K(N). 

gBS _ _j_ ( ^2x2 — 3-2x2 
y/2 V !2x2 12x2- 



Another subset of Sp(2N, R), although not a group, are the transformations of 
the form: 

U(N) = {P\P G Sp(2N, R), P = P T , P > 0} . (2.11) 

These transformations will be called active transformations since they change 
the energy of the system. Generally these transformations are difficult to im- 
plement experimentally, e.g. for states of light strong nonlinear optical media 
are required to realise an active squeezing transformation. 

Example 2 : The one-mode-squeezer changes the variances of 
the quadratures of a state and is an element of H(N). 

gSQ — ( 0\ ^>0is the squeezing parameter 

V d J 

Remark: We will investigate squeezing in the subsequent chapter. 

The only matrix which have both, active and passive transformations, in com- 
mon is the identity matrix. 



2.2.2 Normal Forms of Symplectic Transformations 

We introduce some useful decompositions for symplectic matrices. For a de- 
tailed discussion see Ref. [2] and references therein. 



Polar Decomposition 

From matrix analysis we know that every N x N matrix A can be decomposed 
in two matrices P = V A and U unitary, where for A nonsingular U = P~ l A 
is uniquely determined. If A is a symplectic matrix, then P and hence U are 
real symplectic as well and P £ n(iV) and U £ K(N) by definition. Now for 
all S G Sp(2N, R) exists a unique P £ U(N) and a unique U £ K(N) so that 
S can be written product of both: 

S = PU. (2.12) 

The decomposition shows that every symplectic transformation can be realised 
by a passive transformation, like beam splitters or phase shifters, and an addi- 
tional active squeezing operation. 
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Euler Decomposition (Singular Value Decomposition) 

Since P is real symmetric we may transform it orthogonally to a basis where 
it is diagonal: D = T PO where D and O are again symplectic. For D then 
follows from DaD = a that either a^i = or A^A/ = 1. The diagonal matrix is 
thus of the form D = diag(c? 1 , ...,d,N, g— ) with di > Vi = 1, ...N. 

These diagonal matrices again form a subgroup n(iV) of the symplectic ma- 
trices which we name squeezing transformations or simply squeezers. We will 
introduce squeezed states in Section 3.1 and discuss squeezing as a property of 
covariance matrices in Section 3.3. 

f[(N) = {P\P = diagfdi, d N , — ) with d t > 0} C n(JV) C 5p(2iV, R). 

d± djsi 



We may now write every S G Sp(2N, R) in the non-unique decomposition 

S = U 1 PU 2 (2.13) 

with two passive transformations U\,U2 € K(N) and a squeezing transforma- 
tion P G n(iV). The discrete non-uniqueness results from the possibility to 
arrange the entries of P in different order. 



Iwasawa Decomposition 

The polar decomposition enables us to decompose any symplectic matrix with 
unique factors, but unfortunately one factor is not taken from a subgroup. 
The Euler decomposition allows a decomposition in subgroup factors but is not 
unique. We now introduce a third decomposition which unites both virtues. 



The Iwasawa decomposition states that any S £ Sp(2N, R) can uniquely be 
expressed as a product of three factors 



S = NPU 



(2.14) 



with N G A/", P e U(N) and U G K(N) where M is a subgroup of Sp(2N, R) 
defined as: 



M = { 



A 

C {A~ l ) T 



A 



CA' 1 symmetric > C Sp(2N, R). 



2.3 Moments 

With the postulates of quantum mechanics every system can be described by 
an operator called density matrix or state, usually denoted by p. The density 
matrices provide a more convenient language for quantum mechanical systems 
whose state is not completely known. Additionally every such system can be 
completely characterised by the moments of its probability distribution, as we 
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know it from classical statistics. In the quantum world this is the Wigner func- 
tion, which is in contrast to classical statistics in general not positive. Alterna- 
tively one can use the characteristic function x, being the Fourier transformed 
of the Wigner function and allowing to derive the moments of the probability 
distribution by taking its derivatives. 

Definition 2.2 (Density Matrices) 

The set of density matrices (states) on a Hilbert space 7i fulfilling 

p = fr self adjoint 
p > positive 
tr[p] = 1 normalised 

will be denoted by S(TC). S(TC) is a closed convex subset of the bounded lin- 
ear operators B(7i) on 7i. Note, that all states are by definition trace-class 
operators. 

From now on all p appearing in this thesis are elements of S(7i). The first mo- 
ments of the canonical observables of a state p are collected in the displacement 
vector d and can be calculated via 

d 3 := tr[pR j } = (R j ), j = l,...,2N, (2.15) 

while the second moments of p are collected in the covariance matrix ( CM) 7 
with the relation 

7y := tr[p{Ri - (R z ), R, -{&)}+] (2.16) 
= 2tr[p(R i -(R i ))(R j -(R j ))}-ia ij , i,j = 1, ...,2N. 

The covariance matrix 7 of the state is a real, symmetric, 2N x 2N matrix. 
But not every matrix fulfilling these conditions is a proper covariance matrix, 
from a simple calculation we see that 7 and 7 + ia (and 7 — ia with complex 
conjugation) are necessarily positive. 

Let 7 be the covariance matrix of system of ./V modes. 
Vt; G <D 2N : H( 7 + ia)\v) = ^v*^ + ia^Vj 

= 2tr[pJ2v*(Ri - (^))E^(4- " <4'»] 

* i 
= 2tr[pM Jf M]>0, (2.17) 

with M = Vj (Rj — (Rj)). The right hand side is a positive number for all v 
since p and operators of the form M are always positive. Therefore we find 
that 7 + ia is positive. Similarly one shows that 7 > and 7 — ia > 0. 

It turns out that the restriction 7 + ia > is exactly the uncertainty relation 
for the canonical operators formulated in a basis independent matrix inequal- 
ity (e.g., for all quadratures). Since there are no other requirements than the 
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CCRs and the uncertainty relation, all mentioned properties of 7 are sufficient 
to guarantee that there exists a state p G S(7i) with covariance matrix of 7 . 

Example 3 : The covariance matrix of a one-mode state with 
vanishing displacement is: 

_ / 2(x 2 ) 2(xp) - i 
7 ~\2(px)+i 2(f) 

Assuming (xp) + (px) = 2(xp) — i = 0, the uncertainty relation 
on 7 reads 

/ 2(x 2 ) +i \ ^ n 

Since the trace is positive anyway the only condition to be fulfilled 
is a positive determinant 

4(x 2 )(p 2 ) - 1 > or (x 2 )(p 2 ) > -. 

But this is exactly what we know as Heisenberg uncertainty 
relation for the canonical coordinates. 



Definition 2.3 We define the states of minimal uncertainty to be those states 
where half of the eigenvalues of^A-ia are zero, that is, the uncertainty relations 
hit equality in every mode. 

Example 4 : The condition for the states of minimal uncertainty 
of only one mode can easily be derived : 

7 + i(j=( a . C % I , with eigenvalues Ai + Ao = a + b and 
\ c— i b J 

ab — c 2 — 1 = A1A2. Then for Ai = it is A2 = a + b and thus 
det7 = ab — c 2 = 1 for all one- mode minimal uncertainty states. 
As we will see later (Section 3.1) all states having this property 
are pure displaced squeezed vacua. 



We will need some basic lemmata of matrix analysis to cope with the covariance 
matrices. For an overwhelming variety of such little proofs see Ref. [5] and Ref. 
[6]. Most lemmata used in this thesis are proved in the appendix. 

Lemma 2.1 For every positive N x N -matrix A and every M x N -matrix C: 
CAC^ is again positive. 

Proof : It is W G C M : (v,CACU) = £f m Y% n v* k C kl A lm d n v n = 

Sfm v 'l* ^- lm v 'm — 0) since A is positive. □ 

And finally, how do the moments of a state transform under symplectic trans- 
formations S G Sp(2JV,R)? 
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Lemma 2.2 The covariance matrix 7 and the displacement d transform under 
symplectic transformations S : R 1— > R' = SR according to 

d 1 ^ 5d, (2.18) 
7 ■-> 5 7 S T . (2.19) 

This is easily seen by the definition of the displacement and the covariance 
matrix. The resulting covariance matrix 7' = S 7 S T still fulfills covariance 
matrix properties since SjS T > for 7 > and SjS T + ia = 5(7 + io~)S T > 
for 7 + ia > 0. 

We have done the first steps with covariance matrices which will play an extraor- 
dinary role in this thesis. We will see in Section 3.3 how efficiently properties of 
Gaussian states can be formulated on the level of covariance matrices. Higher 
moments will not appear since Gaussian states are determined by their first 
and second moments, so it is possible to calculate all higher moments out of 
them. 



2.4 Characteristic Function 

The canonical operators x and p have their drawbacks; especially the fact that 
they are unbounded operators on the Hilbert space having "eigenvectors" not 
belonging to it is unpleasant. But there is a possibility to avoid these techni- 
calities. 

Definition 2.4 We define the system of Weyl operators : 

W{ = e^ TaA (2.20) 

for £ E H 2N , a the symplectic matrix and R defined in Section 2.1. 

All are unitary (bounded) operators on the Hilbert space Ti fulfilling the 
following relations: 

Lemma 2.3 (Properties of the Weyl Operators) 

For £,77 S R 2Ar and £ J G 1R 2 belonging to the j—th mode respectively it is: 



tensor product, 



= 1 unitarity, 

W^+rj Weyl relations, 
2N (£, — 77) orthogonality, 

where the CCRs of the canonical operators have been translated to the Weyl 
relations of the Weyl operators. The orthogonality should be taken with a grain 
of salt. 





N 








i=i 


w\ 


= W-c, 




= w\w^ 




= e~^° 
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Proof : All properties are easily calculated except the orthogonality being a bit 
tricky. A physicists proof is given in the appendix, Lemma A.l. □ 

Example 5 : We construct the Weyl operators for a single mode. 
With £ T = (x,p) it follows: W% = exp[i(xp—px)]. In the language 
of ladder operators, defined in Eq. (2.4), we find with 

x + ip „ x — ip 

a = — =— and a = — =— 

the displacement operator D(a) which is often used in quantum 
optics. 

W i{a) = e a * a ~ aat = D(-a). (2.21) 



The name "displacement operator" comes from the fact that the Weyl operators 
generate translations in phase space : 

W^RW\ = ^2e k (e^ TaA R k e-^ TaA ) 

k 

= J~] e k (Rk + €icrij[Rj , Rk] + 0) 

k i,j 

= ^e fc CR fc + &i), 
k 

where we have used the Baker-Campbell-Hausdorff formula, see Lemma A. 2 
and Ref. [7] for a proof. 

The set of Weyl operators is furthermore a basis for the space of bounded linear 
operators B(TC). Therefore it is possible to expand every linear operator from 
B(TL) in that basis with a unique weighting function /a- 

(27t) J]R 2N 

For every trace-class operator A G T l (7i) (e.g., \\A\\\ = tr[\A\] = tr[V At A] < 
oo) we can calculate the weighting function explicitly by using the orthogonality 
of the Weyl operators. We also use that for every bounded B G B{TL) the 
product AB is again in T l {H.), 



tr[AW^} = — f d m r,f A (-r})tr[W v Wz] 

= Ja(0- 

From now on we will use only weighting functions of density matrices which we 
call in analogy to classical probability theory characteristic functions. They 
are connected via: 
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Definition 2.5 The characteristic function of a state p G S(7i) will be denoted 
by 1 Xp an d is given by the expectation values of the Weyl operators 

X P (0=tr{pW ( :}. 

Then the state p of N modes can be expanded in the basis of Weyl operators 
1 



(27T) 



N 



d 2N C x P (-0 w ( . 



R 2JV 



(2.22) 



To be a proper characteristic function, xp has to fulfil the necessary conditions: 



1. 

2. 
3. 



tr[p] = 1 

nt 



x P (o) = i 
p = p' => x P (0 = x* P (-0 

p > => tr[/9AA^] > for all operators A 
especially for an arbitrary combination of Weyl oper- 
ators A = Cfc W^ k with Cfc G C and G IR 2 



,2JV 



for x P follows ^ c fc cf Xp(Cfc - X«) ea^' > 
fe,/ 

4. Xp is continuous at £ = 

For 3. we have used that the trace of a product of two positive operators is posi- 
tive, see Lemma A. 3 for a short proof. The third condition is named a— positive 
definiteness and is the main condition x has to fulfil and surprisingly the above 
conditions are sufficient for x ( see Ref. [8] for a proof). 



The moments of the state p can now be calculated with the help of its charac- 
teristic function: 

{R'kip 



8 2 



£=0 
£=0 



{R'k R i)p 



where ir is the transformed vector of the canonical coordinates i?' = o~R. 

2.5 Unitary and Symplectic Transformations 

We have seen that every state p G S(TC) can be expanded in the basis of the 
Weyl operators = e l ^ TaR with a weighting function x{0 : 

'*-^jL d ™« (2 ' 23) 

Now the question arises how a symplectic transformation on the basis operators 
R change the state pi To answer this question we need an important theorem 
stated by J. von Neumann, Ref. [9], based on the work of M.H. Stone, Ref. 
[10]- 

x We will drop the p if confusions are impossible. 
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Theorem 2.1 (Stone- von Neumann) 

Let and be two Weyl systems over a finite dimensional phase space 

(N < oo), which obey the Weyl relations in Lemma 2.3. 
If the two Weyl systems are 

1. strongly continuous, i.e. G TC : lining \ — W^ll = 0, 

2. irreducible, i.e. V£ G K 2N : [W^ A] = => A oc 1 
there exists a unitary operator U such that 

V£ G R 27V : = u wm r/t. 

When transforming the basis (passive) with a symplectic matrix S: -R i— > R' = 
S R, the state will be expanded in the new Weyl system W sk t 

Our systems of Weyl operators fulfill the necessary conditions for the Stone von 
Neumann theorem, so we know that the two Weyl systems are connected with 
a unitary transformation so that 

p S r = j^y< J^ 2N ^x{-i)u\s)w k ^u{s) 

= U\S)p k U(S). 




Figure 2.1: Passive and active transformations 

When comparing p sk to p k in Eq. (2.23) we see that the above transformation 
could also be done by transforming the state (active) while leaving the basis 
unchanged. In other words to every symplectic matrix S G Sp(2N, R) exists a 
unitary operator U(S) so, that p k i— > p sk = U^(S) p k U(S). 
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2.6 Bipartite Entanglement 

The crucial feature quantum mechanics provides for quantum information the- 
ory is entanglement. It is a purely quantum phenomenon with a lot of counter- 
intuitive consequences. When measuring different parts of a composite system 
(H = Ha Hb) the outcomes of the measurements may be correlated. This 
correlation could be generated by a preparation of the state using only opera- 
tions on each part of the system (local) , classical communication and statistical 
mixing (LOCC). States which could be prepared like this will be called classi- 
cally correlated or separable. But this is not the whole story. In quantum theory 
a stronger kind of correlations can be observed. The so called entangled states 
are non-local in the sense that they can not be prepared by LOCC and that 
both parts of an entangled pair do not have their own properties but contain 
only joint information. 

The definition of separability shows the tensor product structure of locally 
produced states plus mixing. 

Definition 2.6 (Separability) 

A state p is separable with respect to the split A\B iff it can be written (or 
approximated, e.g., in trace norm), with probabilities pi > and YliPi = 
and proper density matrices p^ and p^' belonging to the parties A and B 
respectively, as 

i 

that is, the closed convex hull of product states. Otherwise the state is called 
entangled. The set of separable density matrices p on Ha ®Hb ( with respect 
to the split A\B) will be denoted by S^\g(7i) and is a closed convex subset of 
S(H) defined in Def. 2.2. 

We see that this definition gives rise to complications since it is not easy to 
find out whether a given state p can be written in the form (2.24) or not. A 
significant amount of research has been done to find simpler criteria to decide 
whether a given state is entangled or not. One of the most famous criteria is 
the so called PPT-criterion. 

Theorem 2.2 (PPT criterion for bipartite systems) 

Let p be a state on a Hilbert space H = Ha ®Hb which has a non-positive 
partial transpose then the state is entangled with respect to the split A\B. 

p TA =>■ p is entangled 

In the cases I. Ha = Hb = C 2 , 2. Ha = C 2 and Hb = C 3 and 3. for Gaussian 
states split into 1a x Nb modes this criterion is also sufficient. 

The necessary direction for entanglement was formulated by A. Peres in 1996 
Ref. [11]. The sufficiency was shown by M.+ P.+R. Horodecki (1 and 2) in 
1996 Ref. [12] and R.F. Werner and M.M. Wolf (3) in 2000 Ref. [13]. The 
name PPT criterion is an abbreviation of positive partial transpose. 
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The partial (say A ) transpose of a state is the state transposed only on one of its 
subsystems (only A). In formulae: p a ,n;/3,m = (p TA )p,n;a,m where a, (5 and n, m 
belong to the parties A and B respectively. Note, that p TA is still selfadjoint 
and has trace one. The partial transpose of a state is only a mathematical tool 
and cannot perfectly be performed by a physical transformation. Sometimes 
the partial transpose of a system's state is understood as a time reversal in the 
transposed part only. Phase conjugation in laser beams actually realise such 
a time reversal but it can not be done perfectly since there always comes a 
small amount of noise with it (see Example 7 in the fourth chapter for a little 
discussion) . 

Almost every application in QIT uses non-classical correlations, e.g., telepor- 
tation [14, 15, 16], quantum cryptography [17, 18, 19] and dense coding [20]. 
Because all these funny things strongly depend on the entanglement of the used 
states it is often understood as a resource. Naturally one wants to know how 
much of that resource one has. Due to its importance a great deal of effort has 
been invested in the last decade, to find a sensible measure quantifying entan- 
glement. Several entanglement measures have been proposed, for example the 
entanglement of formation in Ref. [21, 22] and the distillable entanglement in 
Ref. [23] and further discussions in Ref. [24, 25]. We introduce two quantities 
strongly related to entanglement measures, which we will need later on. The 
von Neumann entropy is the the quantum analogue of the Shannon entropy 
Ref. [26] and gives the degree of mixedness or impurity of a state. 

Definition 2.7 (von Neumann entropy) 

The von Neumann entropy of a mixed state p £ S(7i) is defined as 

S(p) := -tr[p lnp] = - ^ Pilnpi 

i 

where the pi are the eigenvalues of the state p. 

The only so far computable entanglement measure is the logarithmic negativ- 
ity, Ref. [27, 28, 29, 30]. With \\A\\ X := tr[y/A*A] we define the logarithmic 
negativity 

Definition 2.8 The logarithmic negativity of a bipartite state p £ is given by 

E N {p) :=ln||p T -*||i. 
with respect to the split A\B. 

If p TA has negative eigenvalues Xi(p TA ) fulfilling 1 = tr[p TA ] = ^ Aj(/? Ta ), 
the sum of the absolute values of the eigenvalues must be greater than one: 
||/9 Ta ||i = j I Aj (p T - A ) I > 1. Hence the logarithmic negativity of p is greater 
than 0. If p has a positive partial transpose p TA > 0, the logarithmic negativity 
vanishes since tr\-\f p Ta p Ta ] = tr[p TA ] = 1. For all separable states this is true. 
But there may also exist entangled states with PPT unless the PPT-criterion 
is also sufficient. In this thesis we will mainly discuss the case where the PPT- 
criterion is also sufficient, namely for Gaussian states of 1a x Nb modes. Note, 
that for all elements of S(H) the logarithmic negativity is positive definite. 
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Chapter 3 

Gaussian States 



In quantum information theory more and more attention is paid to continuous 
variable states. Gaussian states are a nice target to exploit since they are 
available not only theoretically but can be observed and prepared in the lab. 
For example, any laser produces Gaussian states and most optical setup laser 
light can go through preserves this property. From the mathematical point 
of view they are the simplest case of nontrivial CVS showing squeezing and 
entanglement. To learn how one can use continuous variable states for QIT 
purposes Gaussian states play the key role and will be characterised in this 
chapter. 

Definition 3.1 A Gaussian N— mode state p is a state, whose characteristic 
function defined in Def. 2.5 can be written as: 

Xp (0 = exp[-^ T re + iD T C\, e G ^ 2N , (3-1) 

where T is the covariance matrix of the state and D the displacement vector. 

No other parameters appear since in analogy to classical probability distribu- 
tions the quantum Gaussian states are determined by their first and second 
moments alone. When describing Gaussian states, we will often use only their 
covariance matrices since they reflect all important properties Gaussian states 
can have. We will investigate those properties in the subsequent sections after 
introducing some examples of Gaussian states. 

Unfortunately in the literature both, the matrix 7 defined in Eq. (2.16) and 
the matrix T appearing in the characteristic function above, are called co- 
variance matrix. Usually it is clear which one is meant and therefore we will 
not distinguish them explicitly. As a general rule in mathematical formulae we 
will use small greek letters in the first case and capital letters in the second. 
The same holds for the displacements d and D. The two displacement vectors 
and the two matrices can be transformed into each other using the symplectic 
matrix defined in Eq. (2.2) with the following transformation law: 

D = ad and T = a 7 a T . (3.2) 

As one can easily check the T-matrix again fulfills covariance matrix properties. 



19 



20 



CHAPTER 3. GAUSSIAN STATES 



3.1 Coherent, Squeezed and Thermal States 

We introduce three important classes of pure one mode Gaussian states, namely 
the coherent, the squeezed and the thermal states of an electromagnetic field 
mode. General reference for this section is Ref. [7] by S.M. Barnett and P.M. 
Radmore. 



Coherent States 

The coherent states are defined as the eigenvectors of the non hermitian anni- 
hilation operator a : a\a) = a\a), with a being a complex eigenvalue, and one 
easily finds that those state vectors are superpositions of number state vectors 
of the form 

They can be generated from the vacuum with the unitary operators D(a) in- 
troduced in Eq. (2.21) which induce displacements in phase space. 



D(a) = e a ^- a * & = e a ^e~ a * & e-2, 
\a) := D(a)\0). 

They are an overcomplete nonorthogonal set of state vectors spanning the entire 
Hilbert space resolving unity, 

1= — —\a)(a\, 

and the scalar product between two coherent states 

(a\(3) = 

With a = r?1 ^| r?2 and the identification Eq. (2.21) we calculate the characteristic 
function of the coherent states: 

Xa (0 = tr[\a){a\Ws]=tr[D(a)\0)(0\D(-a)Ws] 

= tr[W-r,\0)(0\W v W^ = e~^ T < ^W^W-^) 

= e -f" T ^ e 5(^) T -'? (0|W^|0) = e- ir > T °S 

= e -iv^(»(a),3(a))a{ 



where we used the cyclicity of the trace, the Weyl relations from Lemma 2.3 
and (0|W^|0> = e'^^. 



The coherent states are thus of Gaussian type and furthermore all of them 
have the same covariance matrix: T a = 1, while the displacements depend on 
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the value of a: D r 



V2 



9(a) 
-3? (a) 



Similarly, we have j a = aT a a T = 1 



and d a 



-aD a = V2 



9(a) 



. The uncertainty relation is exactly fulfilled 



with A 2 x • A 2 p = j and is equally stretched in all directions in phase space. 
Therefore coherent states can be pictured as circles of diameter -k= around the 

points (v25R(a), \/29(a)). For a = the coherent state vector |0) is just the 
vacuum with Tq = 70 = 1 and Dq = do = 0. 




Figure 3.1: States in phase space 

The picture shows the vacuum state (blue), a coherent state (yellow) and a 
displaced squeezed vacuum (red). They all span the same area meaning that 
the product A 2 x ■ A 2 p is the same for all states, namely it is \, the lower bound 
of the uncertainty relation. Hence all these states are minimal uncertainty 
states. 



Squeezed States 

Similarly one can produce states with non-energy conserving squeezing opera- 
tors S'(C)) though experimentally it is difficult to realise such a transformation 
requiring strong Kerr media. 

S(C) := e -§ a2t + ^ 2 

— 77 cx.p[i(f)] tanh r a?^ — \ lnfcoshrVa^a+da'h \ exp[— tanh r a 2 

S is unitary since it is 5(C) = 5t(-C) with C = re**. The squeezed vacuum is 
prepared by letting a squeezing operator act on the vacuum. 

10 := £(010) 



v 7 cosh \ 



(~ ex.p[i<f>] tanhr)" |2n). 



n=0 
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The characteristic function of the squeezed vacuum with a = ^~^ 2 is: 

= (0\S(-C)D(-a)S(O\0) 
= {0\S(-()S(()D(-a')\0) = e~'4- 
with a 1 determined by commutator of D(—a) with «S(C) and given by 
, £1 [cosh r + cos (J) sinh r] + £2 sin <f> sinh r 

a = 

£1 sin <j> sinh r + £2 [cosh r — cos sinh r] 
+ * 71 

we get 

t^(£) g — (wsh 2r+sinh2r cos (^)) g—j^| (cosh 2r— sinh 2r cos <^i) 

x e - ^ 1 ^ 2 shlh2rsin ^ 

and hence all squeezed vacua are of Gaussian type. They have vanishing dis- 
placements and covariance matrices 



cosh 2r + sinh 2r cos 4> sinh 2r sin <f> 

sinh 2r sin rf> cosh 2r — sinh 2r cos 1 



From the CM one can read off that the variance of the state in one direction of 
the phase space is different than in the other; that is why the states are called 
squeezed. But the uncertainty relation is still fulfilled as a short check shows. 

A 2 x ■ A 2 p = — (cosh 2r + sinh 2r cos (j)) (cosh 2r — sinh 2r cos 4>) 

1 2 o 1 sinh 2 2r sin 2 6 1 

= —(cosh 2r — sinh 2rcos<p) = — I > — , 

4 4 4 4 

where equality (minimal uncertainty) is reached for r = or eft = 0. In the first 
case there is no squeezing at all, S(0) = 1, and we again have the vacuum. In 
the second case has the form 

/ cosh 2r + sinh 2r \ _ / e 2r 

C=r ~ ^ cosh 2r - sinh 2r )~\ e" 2r 

where we recognise that r^ =r could be prepared from the vacuum with the 
symplectic squeezing transformation introduced in Example 2 of the second 
chapter with the squeezing parameter d = e r . All squeezed states can be 
prepared from the vacuum by a squeezing operation (squeezed vacuum) and an 
additional displacement in phase space. 
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Thermal States 

The statistical concept of equilibrium states or Gibbs states is naturally taken 
to the quantum world. We define the thermal state of a system described by a 
Hamiltonian H as the exponential of the Hamiltonian together with the factor 
(5 = 75— y, where ks is the Boltzmann constant and T is the temperature of 
the system. For a one mode harmonic oscillator the well known Hamiltonian is 
given by H = u>(n + with h = aJa being the number operator of the mode 
and h = 1 throughout this thesis. The thermal state for this Hamiltonian is 
given by 



Pt 



-(3H 



tr[e-P H ] 



£ 

n,m 



n){n\ 



-/3u;(n+i) 



\m){m\ 



trie 



-/Mn+|)l 



n)(n\, 



with tr[e-^ n+ ^] 



2sinh^ 



The characteristic function then is: 



^|n)(n|e-^ n+ 5)# € 



Xt(0 = 2 sinh^ tr 
Xt (0 = 2sinh^ ^e-^ + D(n|^|n). 

n 



The elements of the Weyl operators in the number basis with a 



6 +''€2 
V2 



are 



(n\W^\m) 



(n\D(-a)\m) = (0\-^D(-a ' 



f m 



10) 



(-H 2 )' 



I a I 

^ (-a)"" m v 7 ^! ELo n { m-m+n- m y. 
ya ) ym-ii- 2^i=q i\(n-i)\{i +m - n )\ 



(-a)»- 



^L^- m (|a|^) m<n 



; 2 (a*) 



m— n / re! Ttn—n 
m! ^re 



l (|a| 2 ) m >n, 



m < n 
m > n 



where L^(x) are the Laguerre polynomials; so we get for n = m : 



(n\Wz\n) = (n\D(-a)\n) = e"^"L° (|a| 2 ), 
and Xt becomes 

Xt(0 = 2sinh^e-^^e-^^)L0(|a| 2 ). 
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With sum formula (8.975.1) in Ref. [31] we finally get 
Xt(0 = 2smh— e 2 e 2 - — — e . 



2 e 



1 



*l+*2 



Hence the thermal state is a Gaussian state with vanishing displacement and 
covariance matrix 

1 



tanh 4jr 



1 = 7t, (3-4) 



that is, a diagonal 2x2 matrix with identical entries. We will meet the thermal 
states again when discussing normal forms of CMs in the next chapter. In the 
limit (3 = -^Fp — ► 00 we get tanh 4^ — > 1 and hence 



X(T=0)(O = 

that is, the thermal state becomes the vacuum in the limit of zero temperature. 



3.2 Normal Forms of Covariance Matrices 

Covariance matrices are real symmetric matrices, hence it is always possible to 
diagonalise a CM with orthogonal matrices to the standard normal form which 
is the diagonal matrix of the eigenvalues of the CM. But other, symplectic, 
normal forms exist which are more convenient for covariance matrices because 
the resulting quantities are invariant under subgroups of the symplectic trans- 
formations. We introduce two very important normal forms which are used 
extensively in the following chapters. 



3.2.1 Williamson Normal Form 

Surprisingly it is possible to diagonalise real positive symmetric and even- 
dimensional matrices with symplectic matrices. The resulting diagonal matrix 
is called Williamson normal form and the proof of existence and uniqueness 
was given by J.Williamson in 1936, see Ref. [32]. 



Theorem 3.1 (Williamson normal form) 

For every real strictly positive symmetric 2N x 2N matrix A there exists a 
matrix S £ Sp(2N,H), so that 



A WNF = SAS T 



( ai 

01 



\ 



\ 



ON 

a/y ) 



(3.5) 
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with the symplectic eigenvalues en > for all i = 1,...,N. The symplectic 
matrix S and the diagonal entries ai are unique up to permutations of the 
diagonal entries in A WNF . 



Remark: If A had a vanishing eigenvalue it would not be possible to bring 
it to the above form. Instead, Williamson-diagonalise the regular, positive 
principal submatrices of the matrix A. Now transform the singular subma- 

trix Aq of A to its orthogonal diagonal form Dq = O t AqO = ( n ] with 

V a o J 

O G 50(2, R) C 5p(2, R). The symplectic eigenvalue of Aq is counted as zero. 

Note, that the determinant of A is easily calculated from its symplectic eigen- 
values: det A = J\i=i a i an d that the symplectic eigenvalues of a matrix are 
invariant under symplectic transformations. 

Especially covariance matrices can be Williamson-diagonalised and as we will 
see in the following discussion of the Williamson normal form the transformed 
covariance matrix Y WNF = ®^ 1 s«l2x2 describes the state of N uncoupled 
harmonic oscillators in a heat bath with temperatures Tj = Tj(sj). 



Symplectic Eigenvalues 



In general it is difficult to find the symplectic transformation which brings a 
given matrix to Williamson normal form, but the eigenvalues are easily calcu- 
lated with the following lemma 

Lemma 3.1 The symplectic eigenvalues of a real symmetric positive 2N x 2N 
matrix A are the positive eigenvalues of the matrix ioA. 



Proof : With A = S- 1 A WNF S^ T 
spec(iaA) 



spec(fa5- 1 ^ H/iVF 5- T ) = spec(iS T aA WNF S~ T ) 
S pec(iS- T S T aA WNF ) = spec(iaA WNF ) 

N 



spec i 



3=1 



1 
-1 



were we used the properties of symplectic matrices and the cyclicity of the 
spectrum, see Lemma A. 4. The eigenvalues Aj of iaA WNF are calculated via 



= det[iaA 



WNF 



Ail] 



N 

n 

3=1 



-Aj iaj 
id* — A,- 



N 

© 

3=1 
N 

3=1 



-A,- 



-zaj 



we find the spectrum of io~A to be spec(io~A) = {a{, —a{}^ =1 completing the 
proof of the lemma. □ 
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3 1 \ 

Example 6 : Given a covariance matrix V = I I then 



1 1 / 

the symplectic eigenvalues of T are the positive eigenvalues 

of iaT = ( \. 1 . ) . The eigenvalues are Ai 2 = ±\/2, so 
y — 3i —1 ) 

the symplectic eigenvalue of V is T± = y2. Even faster is 
detr = 2 = rf, thus Ti = sft. 



Thermal States 

Question: Given a covariance matrix in Williamson normal form and displace- 
ment D, what is the corresponding Gaussian state p ? 

We first see how block diagonal covariance matrices T = ©j = i Tj are related to 
product states. 

Lemma 3.2 (Block diagonal CMs) 

A Gaussian state p with a covariance matrix in block diagonal form is a product 
state, and conversely if p is a state of product form then its covariance matrix 
is of block diagonal form. 

Proof : We use the decomposition of p in the Weyl system according to Eq. 
(2.22) and insert the characteristic function with the given block diagonal ma- 
trix T = 0^ Ti and displacement D. For £ G R 2N and & G R 2 , Vi = 1, ...N, 
p can be decomposed 

N 

= ® 

i=l 

iV 

i=l 

where the /9j are one-mode Gaussian states with covariance matrices Tj and 
displacements -D,,. □ 

We learned in Section 3.1 that the one-mode thermal states of the harmonic 
oscillator have the covariance matrix Tj = - - (iuj 1 and vanishing displacement. 

Since there is a one to one correspondence between the pairs T and D and 
the Gaussian states p, the displaced thermal states are the only ones having 
such a diagonal matrix. The Gaussian state belonging to a CM in Williamson 
normal form F = Sjl with symplectic eigenvalues Sj and displacement D 

is a product state of displaced thermal states with independent heat baths of 



J d%e^^W^ 
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temperature Tj. 



A' 



N 



P 



Pi 



W nP t{Ti)W 



i=l i=l 

The symplectic eigenvalues of the CM are then connected to the temperature 
of the heat bath namely 
1 



tanh 



2k B Ti 



3.2.2 Simon Normal Form 

We introduce a normal form, as proposed in Ref. [33], acting only locally on 
the first and second mode respectively and therefore preserve the entanglement 
properties of a given state p. 

Theorem 3.2 (Simon normal form) 

For every two-mode covariance matrix V there exist matrices S\ G Sp(2, R) and 
S 2 £ Sp(2, R) so that: 

fa 

(Si®s 2 )r(s[ @sr 



^SNF 



C 

V o 





a 

d 



d 


b ) 



This normal form is unique up to the little subtlety that only the relative sign 
of c and d is determined. 

The transformations S = S\ ® S 2 are local transformations acting only on 
the first and second mode respectively, in particular they do not change the 
entanglement properties of the covariance matrix they are acting on. 



Simon Invariants 

The determinants of the submatrices of V and V itself are preserved by the 
diagonalising symplectic transformation and appear as the Simon invariants 
a, b, c, d with the relations: 

A C 
C T B 



j a = V det A, 



b = VdetB, (3 - 6) 
detr = a 2 b 2 + c 2 d 2 - ab(c 2 + d 2 ) cd = detC. 

We hence reduce our ten parameter covariance matrix to four real numbers still 
reflecting the entanglement properties of the original CM. 



3.3 Properties of Characteristic Numbers 

All derived properties of the covariance matrix 7 translate identically to the 
matrix T = a ^7 a 1 \ Namely T is a real, symmetric, positive matrix fulfilling the 
uncertainty relation and has the same determinant, eigenvalues and symplectic 
eigenvalues as 7. 
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Simon Invariants and Symplectic Eigenvalues 

In the last section we saw that there exist normal forms respecting the symplec- 
tic character of the canonical transformations on covariance matrices. These 
forms were related to invariant quantities whose properties we will discuss in 
this section. We will also see that properties of the states like squeezing and en- 
tanglement are strongly connected to these numbers. We denote the eigenvalues 
of a given N— mode CM V by Aj, for i = 1, ...,2N, the symplectic eigenvalues 
by Tj, for j = 1,...,N and the Simon invariants (only for two- mode systems) 
by a, b, c, d. 

Assume we have a covariance matrix V in Simon normal form and we want to 
know the eigenvalues and the symplectic eigenvalues of Y SNF '. We find that the 
eigenvalues of T SNF , expressed with the Simon numbers, are given by 

2 2 (3.7) 

x SNF = a_ + b ± l_^ {a _ h)2+M2 

Remark: These numbers in general are not the eigenvalues of V since in gen- 
eral the symplectic transformations Si and S2 are active transformations (see 
Section 2.2). 

Furthermore, the symplectic eigenvalues of a two mode system with covariance 
matrix Y SNF can be expressed by the Simon invariants via 



1 



r lj2 = — \J a? + b 2 + 2cd ± v^a 2 + b 2 + 2cd) 2 -4det 7 . (3.8) 
V ^ 

and are the same as the ones of T itself. 

Positivity From the positivity of T it follows that 
Xi > 0, i = l,...,2JV, 

T 4 > 0, j = l,...,N, (3.9) 
a, b > 1 and c 2 < ab, d 2 < ab. 

Uncertainty Relation The Heisenberg uncertainty principle formulated in 
Section 2.3 is also true for the capital matrices T (see Eq. (3.2) and Lemma 2.1). 
Since the matrix T + ia is positive it can be brought to Williamson normal form 
by symplectic transformations S without changing the symplectic eigenvalues. 
We have: 

S [T+ia] S T = STS T +i a = T WNF +ia = ( V J>. ^ \ . (3.10) 

j=l V 1 3 J 

To get a positive matrix all eigenvalues /ij, i = 1, ...,2N, of the above matrix 
have to be positive. The equations determining these eigenvalues are 

N 

11 [(Tj - fa) 2 - 1] = 0, for % = 1, 2N, 

3=1 
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and hence 

Tj ± 1 = m > 0, Vi = 1, iV, Vj = 1, 2iV 

Thus the restrictions on the symplectic eigenvalues of T are 

Tj > 1, for all j = 1,... ,iV. (3.11) 

Note, that the determinant of T is therefore always greater than or equal to 
one. For the Simon invariants we find the following, only necessary, criterion 

detr>det^4 + det-B + 2detC + l, (3.12) 

where A,B,C are the submatrices of V connected to the Simon invariants via 
Eq. (3.6) 

Pure States 

Definition 3.2 A state p is called pure iff it can be written as a projector 
p = \ifj)(ijj\ whereas a mixed state is a convex combinations p = Pi\ipi) {"4>i\ 
of the states \ipi){ipi\ with probabilities pi (pi > 0, J2iPi = ^) an ^ a t least two 
i so that pi > 0. 

It is well known that a state p is pure iff tr[p 2 ] = 1, but how can one express 
purity in covariance matrix language? 

Lemma 3.3 A Gaussian state p with CM T is pure if and only if 

detr = 1 : Fi = 1 (ioT) 2 = 1, 

where the Tj are the symplectic eigenvalues of T. In the two-mode case the 
conditions on the Simon invariants are 

a = b = coshr and c = —d = sinhr, 

where r is a squeezing parameter of the CM. 

Proof : We expand the state in the Weyl operator basis and calculate the trace 
of its square 

tr[p 2 ] = J^ym J I ^riXpi-OXpi-^tr^W,] 
1 

(2^p 
1 

With p being a Gaussian state with covariance matrix V and arbitrary displace- 
ment we get 



d 2N cJd 2N r, Xp (-0x P (-v)S 2N (C + v) 

d 2N t x P (-0x P (0- 
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Hence p is pure iff detT = 1. It follows directly that the symplectic eigenvalues 
Ti, all necessarily greater than or equal to one, just take their lower bound. 
Furthermore half of the eigenvalues of iaT are one while the others are minus 
one. Hence (ioT) 2 has all eigenvalues equal to one and is the identity matrix. 
For the Simon invariants we easily get the above relations by inserting the Tj 
in Eq. (3.8). □ 

Squeezed States If a (general) CV state has a smaller covariance in one 
direction of the phase space than in another it is said to be squeezed. For a 
basis independent formulation one has to consider all those symplectic trans- 
formations which do not change the degree of squeezing. We encountered those 
transformations when discussing normal forms of symplectic matrices in Eq. 
(2.10), where we learned that they are in the intersection of the symplectic and 
the special orthogonal group. 

Definition 3.3 An N—mode state p with covariance matrix T is squeezed if 
and only if 

3S e K(N), 3k : (STS T ) kk < 1. 

Here K(N) = Sp(2N,H) H SO(2N, R) are the passive symplectic transforma- 
tions. 

We picked only those states which have a smaller covariance in at least one 
phase space direction than the vacuum state. Unfortunately, the criterion is 
not practicable, but it can be formulated in terms of the eigenvalues of the CM 
as the following lemma shows: 

Lemma 3.4 An N—mode state p is squeezed if and only if the smallest eigen- 
value \ m in of the state 's covariance matrix is smaller than one. 

Proof : First we show that the diagonal entries of a symmetric matrix A 6 
Mat (AT x 2V, R) are always bigg er or equal than the smallest eigenvalue arnin 
of the matrix. Let {a k } k= i be the orthonormal system of eigenvectors of 
A and {ej}jLi an arbitrary orthonormal basis in H N , then the diagonal en- 
tries of A in that basis are An = (ei\A\ei) = Y,k,i=i (zi\ a k) (ak\A\a{) (ai\ei) = 

Ya=i a i\( a i\ e i)\ 2 > a min YaLi \( a i\ e i)\ 2 = a min, as stated. The passive trans- 
formations in K(N) do not change the eigenvalues of the CM and since 
all diagonal entries are bigger or equal to the smallest eigenvalue we have 
VS £ K(N), V/c : (S~fS T ) k k > ^min with equality when S diagonalises V 
to its orthogonal normal form and the k— th entry is \ m in- Hence the state can 
only be squeezed if and only if the smallest eigenvalue of the CM is smaller 
than one. □ 

Entangled States In Section 2.6 we defined what entanglement is and how 
one can decide whether a given state p is entangled or not. We introduced 
the partial transpose of a bipartite system and stated the PPT criterion for 
states p £ S(Tl). For Gaussian states it is possible to formulate this criterion 
on covariance matrix level. We first see what effect the partial transposition of 
a state has on the states covariance matrix. 
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Lemma 3.5 The partially transposed covariance matrix (on part A) of a bi- 
partite state p of N = N A + Nb modes with covariance matrix T is defined by 
the state p TA with the covariance matrix T Ta , which is given by 

T Ta = (-M A et B )T(-M A et B ), (3.13) 



with M A = M T A = ® Na 



i=l 



1 

-1 



Proof : Partial transposition is a basis dependent non-unitary operation and 
since it is not a completely positive map, it can not be realised perfectly by 
physical means. Mathematically the transformation can be done, since the 
spectrum of the partial transpose is not basis dependent, but we have to choose 
a basis for the calculation. We denote the (multimode) number basis vectors 
by \n A ) and |n#) for 7i A and TCb respectively and calculate the characteristic 
function of the passively transformed state p i— > p Tjv : 

X P (0 = tr\pWfa] = ^2 ( n AnB\pW^\n A n B ) 

n A ,n B 

= ^ { n AnB\p\m A m B )(m A \Wf {A ^ A \n A )(mB\Wf iBiB \n B ) 

n A ,n B ,m A ,m B 

= ^2 ( n AnB\ p\m A mB}(n A \W^Jm A }(m B \Wf i , g ^ B \nB} 

mA,n B ,nA,m B 

= tr\p TA W k ,^. 



To determine the passive transformation on the basis operators R, bringing p 
to p TA , the above expressions shall be equal. We assume that on system B 
no changes take place, Rb = R B i while for system A a linear homogenous 
transformation 

N A 

M A = M fulfilling R' A = M A R A 

i=i 

has to be determined. For every single mode i € A we have by the definition of 
the Weyl operators: 

{milWfaJm) = {n i \W M ^m i ) = {n i \e i ^ aMki \rn i }. (3.14) 

The 2 x 2-matrix M is not a symplectic matrix because if it was one, partial 
transposition could in principle be realised perfectly. Additionally from sym- 
metry considerations M has to be selfinverse so that the determinant of M 
has to be minus one. Therefore it fulfills MaM T = —a and expression (3.14) 
becomes: 

(ni| e ^ T<jM ^|mi) = He-^'^^K) 
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With the complex numbers cti = and 0i = Yl"j=i lj ^7= 2j - and with 

Eq. (3.4) we get : 

(mi\W^.\ni) = (mi\D(-a)\ni) 

e-^(- ai ) m >- n >^L™*- n i(\ ai \ 2 ) ni < m t , 
e-^(a*) n >- m *^L%- m >(\ ai \ 2 ) m > m u 



\0i\ 2 



(mi\W AiM ^,\n 



hence = a* and therefore M 1 = M = ^ ^ ^1^) = ^ T ' ^ n ^ e canonical 

operators, the displacement and the covariance matrices, this transformation 
acts according to Eq. (2.7) and Eq. (2.18) (even if the transformation is not 
symplectic but linear homogeneous): 

R^R' = (Ma ©1b) A, 

7^7 Ta = (Ma © Is) 7 (Ma © 1b) > (3.15) 

r^r TA = (-M4©i B )r(-M4ffii B ), 

since ctM^ = — Mac □ 

Note, that partial transposition leads to a reversal of all momenta in the system 
A while the positions stay unchanged and the system B is untouched. We finally 
got a recipe of calculating the partially transposed covariance matrix. 
Attention: There is a little inconsistency with the definiton of the partially 
transposed covariance matrix since one does not transpose one part of the ma- 
trix, as one does when calculating the partial transpose of a state. Instead, it 
is the covariance matrix associated to the partial transpose of the state and is 
calculated by appling a non-symplectic transformation T = — Ma ©1b- 

We easily see that Y T - A is still a real, symmetric, positive matrix but if the partial 
transpose of p fails to be positive its covariance matrix T Ta by Eq. (2.17) does 
not fulfil the uncertainty relation anymore and vice versa. Finally the PPT 
criterion on covariance matrices is formulated in the following theorem. 

Theorem 3.3 (PPT criterion for bipartite Gaussian systems) 

Let p be a state on a Hilbert space Ti = TCa © with a partially transposed 
covariance matrix T Ta not fulfilling the uncertainty relation, then the state is 
entangled with respect to the split A\B. 

T Ta + io~ =>■ p is entangled. 



In case p is a Gaussian 1a x Ns-mode state this criterion is also sufficient [13]. 
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As we saw in Eq. (3.10) the uncertainty relation for a matrix is equivalent 
to restricting the symplectic eigenvalues of the matrix on values greater than 
or equal to one. We summarise that any state with covariance matrix T is 
entangled if its partially transposed matrix T Ta has an eigenvalue smaller than 
one. 

T TA +ia^0 <^ 3 j- T J A < 1, ( 3 - 16 ) 

with tJ a the symplectic eigenvalues of F Ta . It follows immediately that at least 
one of the (usual) eigenvalues of T Ta is smaller than one and since spec(r TA ) = 
spec(r) every entangled state is necessarily squeezed. 



Lemma 3.6 The logarithmic negativity, defined in Def. 2.8, of a Gaussian 
state p and CM T can be calculated with the help of its partially transposed 
covariance matrix T Ta having the symplectic eigenvalues Tj A : 

Em(p) = - E ruin (hidf 4 ], 0) =: E X (F) =: E^j). (3.17) 

i 

Proof : As we saw in Section 3.2, every Gaussian state can be brought to a prod- 
uct of thermal states with different temperatures by applying an appropriate 
displacement and symplectic transformation. The state was then written as 



oo 

£ 



Pm \ni){rii 



where the probabilities are given by p ni = 2sinh e ft a; (™!+2) which are con- 



nected to the symplectic eigenvalues Tj of the CM of p via Tj 



tanh ■ 



or after 



a little calculation p ni 



The construction also applies for the 



Fi+i \r i +i / 

partial transpose of a state but since p Tjv is not necessarily positive the prefac- 
tors in p^ A are not probabilities anymore, since they can be negative although 
they still sum up to one. The relation 



P 



T A 



1 



rJ A + i 



rf A + i 



still holds, but it is not possible to assign a sensible temperature to the state 
p T A _ "Yy e calculate the logarithmic negativity of a Gaussian state p: 



E N {p) = In 1 1/- 4 1|! =ln 




Pnt \ni)(ni 



where we used that the diagonalising unitary transformation on p Tjv leaves the 
trace norm invariant. 
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Applying the definition of the trace norm gives: 



^2(Pn?) 2 \ni){ni\ 
■i \m=o 



In tr 



Yl \Pnf \ \ n i){ n i 



i \n;=0 

and inserting the p^ A gives 



n 



E 1*1 

.m=o 



E N {p) = £ fin 



.n^o rf* + 1 



rT A - i 



rf 4 + i 



rf- 4 + i - ir;- 4 - i 



- -e{^, |^=-e^(^].o). 

The same logarithmic negativity results when using 7 Ta , having the same sym- 
plectic eigenvalues as T Ta . □ 



3.4 States of Maximal Entropy 

Definition 3.4 T/ie logarithm of a state p £ S(7i) is defined via its diago- 
nalised normal form diag (pi,p2, ■■■) = U' pU according to 

ln(p) := U diag(lnpi, lnp2 ...)W 

where p±,p2,--- are the positive eigenvalues of p and lnpj the logarithms of the 
eigenvalues. 

Lemma 3.7 The logarithm of a tensor product of two states pa and ps is given 
by 

ln(p A ® p B ) = ln(pA) ® 1b + 1a ® m (Pfi) 
Proof : The tensor product of two states can be diagonalised locally with 
(U A ® U B ) (pa ® pb) [U\ ® U B ) = diag(pf ,p£, ...) (g> diag(pf ,pf , ...) 



A R 
PlP2 



A R 
P2Pl 



P2P2 



dim 

[pfdiag(pf,pf,...)] 
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where pf and p B , are the of eigenvalues of pA and ps respectively. The loga- 
rithm of the tensor product is by Def. 3.4 given by 

/ lnpf + lnpf 



ln(p A (g) p B ) 



{Ua ® U B ) 



lnpf + lnpf 



B 



lnpf + lnpf 

lnpf + lnpf 



V 

/ lnpf i s 

Inputs 



(4® 4) 



V 

/ lnp^ 



B 



+ (U A ® 17b) 




lnpf 



lnpj; 



B 



lnp 



V 



(4 ® 4) 



C/a diag(ln pf, lnpf, ...) (8) is 



+ii ® 17 B diag(lnpf ,lnpf /s , ...) 17+ 
= ln(p A ) ® 1b + 1U ® ln(ps) 
as stated in the lemma. 



□ 



Theorem 3.4 (States of Maximal Entropy) 

Of all states with fixed first and second moments, the Gaussian state maximise 
the von-Neumann entropy. 

Proof : The definition of the von Neumann entropy was given in Def. 2.7. We 
again transform the Gaussian state p to a product of thermal states 

N 

p< = u ] P u = 

1 = 



where 77 



U3 



. The logarithm of p' is with Lemma 3.7 given by 



N 



lnp' = J2h 



-1 <S> 



-/3iUJ 



In (tr[e- PiAi ] \ 1 



1n- 



where all addends shall be understood as identity operators on all N — 1 modes 
except the i-th mode where the above operators are inserted. For the Gaussian 
state p and an arbitrary state 5 with the same first and second moments 77 and 
r we calculate the difference of their entropies 

S(p)-S(6) = S(6\\p)+tr[(5-p)lnp] 

> + tr[U(5 - p)UUnp] 

= tr[(5' — p) In p'\. 
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Here 5' and p' are the transformed states with vanishing displacement and 
covariance matrix Y WNF and S(p'\\p) the relative entropy of p and p' which is 
a positive quantity, see proof in Ref. [34]. Finally tr[(5' — p') In p'] = since from 
the above discussion we see that In p' is a polynomial of second degree in the 
canonical operators i,i and pi (i = 1, . . . , N), and hence picks out only the first 
and second moments of the difference, and because 5' and p' have identical first 
and second moments, they vanish. This completes the proof of the theorem. 
□ 



Chapter 4 

Gaussian Operations 



4.1 General Gaussian Operations 

Quantum operations preserving the Gaussian character of all Gaussian states 
are naturally called Gaussian operations, see Ref. [35, 36, 37] and Ref. [38, 39, 
40] for Gaussian channels. As we want to stay in the set of Gaussian states we 
have to take care that the operations we apply do not drive us out of that set. 
In the following we will briefly discuss the classes of Gaussian operations. 

Gaussian Unitary Operations The unitary evolutions U = e lHt with the 
Hamiltonian H which preserve the Gaussian character of all Gaussian states 
are those where H is a polynomial of second degree in the canonical operators 
Xi and pi for % = 1, . . . , N. One understands that by recalling that a Gaussian 
state p can be brought to an exponential form were the exponent is a polyno- 
mial of second degree in the basis operators (see the proof of Theorem 3.4). 
Only a Hamiltonian which is again such a polynomial can transform all Gaus- 
sian states to Gaussian states. We immediately see that these Hamiltonians 
can only generate translations in phase space and symplectic transformations 
S £ Sp(2N, R). The most relevant operations which we can implement experi- 
mentally are of Gaussian type. An example of a non-Gaussian unitary operation 
is the realisation of the Kerr effect whose Hamiltonian is proportional to third 
powers of the ladder operators. 

Gaussian Dilation / Channels Consider an N — mode Gaussian system cou- 
pled to the environment it is living in. In covariance matrix language the CMs 
of the system and the environment sum to the covariance matrix of the whole, 
according to 

7kj = 7s ©7e- 

When applying a transformation on the system we have to take into account 
that the system and the environment always interact with each other and the 
environment may evolve during this process. We assume that this interaction 
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is of symplectic type. Hence the CM of the whole gets transformed by 

7™ i— > S j w S T with S = 

with S s (S e ) the part of the transformation belonging to the system (environ- 
ment) only and the Si describing the interaction between the system and the 
environment. Note, that the submatrices of S are not necessarily symplectic. 
Since we only observe the system the behaviour of the environment is neglected. 
In mathematical formulae we take the average over all possible configurations 
of the environment, that is on states, we trace out the environment. 

p' s = tr e [U{S) Pw U\S)]. 

with U(S) being the unitary transformation associated to the symplectic trans- 
formation S. The corresponding mathematical operation for the covariance 
matrices of the states is to take the principal submatrix of the covariance ma- 
trix of the whole, belonging only to the system. This reduction is denoted by 
brackets with an index [ ] s . 

is = [s-y w s T ] s . 

The resulting covariance matrix for the system alone is then 

% " 'U, 2 S e ){ je)\Sf tl St )\ s 

= SsJsSj + Si^^feSf^. 

With the coupling we do not get only the desired transformation S s 7 S Sj but 
also an additional part depending on the state the environment is in and the 
interaction between system and environment. Formally the symplectic trans- 
formation can be done on the composite of the system and the environment. In 
reality that happens automatically and with our ignorance of the environments 
degrees of freedom we only can see the effects on the system, mathematically 
formulated by taking the reduction on the system while neglecting the envi- 
ronment. We recognise that the interaction between system and environment 
introduces some noise in the system denoted by G = S^i 7 e Sf 1 . Thus, the 
coupling leads to a decoherence process in the system. 

But can we not prepare everything so that the noise vanishes? We have to make 
sure that S is a symplectic matrix and that the resulting 7 S has still covariance 
matrix properties. So first of all G has to be symmetric what it absolutely is. 
Secondly we have S s cr s S s + Si : i a e Sf l = a s and thirdly the uncertainty relation 
has to be fulfilled by 7 e : 

7e + i&e > => 7 e Sf A + % G e Sf^ > 

G = S it i 7 e 5^ > — i S it i a e Sj^ 
G + ia s — iS s CF s Sj > 0. 

We see that the noise G has to fulfil a kind of uncertainty relation as well, 
depending on the transformation S s we want to implement on the system. 



Ss 
Si,2 



Si,i 

Se 
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Especially the case G = is allowed only in case S s was itself a symplectic 
transformation. We conclude that all experimentally available operations, not 
only those of symplectic type for the system, can be done but there is always a 
quantum lower bound for the precision of the transformation which can not be 
beaten. 

Example 7 : We determine the minimum noise for a time rever- 
sal, e.g., phase conjugation of laser light. The transformation on 
the covariance matrix of the system is done with Ma , the matrix 
for partial transposition or time reversal introduced in Lemma 3.5 



7 i-> M A 7 M A = 




Thus the noise G has to fulfil the uncertainty relation 

° <- -ok; -°o 

< G + 2ia. 

As we learned in Section 3.3 such an uncertainty relation requires 
symplectic eigenvalues of G greater or equal to two. A perfect 
phase conjugation is thus not possible, the precision can be 
arbitrary small when using, e.g., strong laser pulses, but we can 
not reduce the noise completely since we have to fulfil a minimal 
quantum limit. 

Adding of Classical Noise Adding classical noise to a Gaussian state p 
is described by a convex combination of random displacements of the state p, 
distributed with a normalised Gaussian distribution 

with a positive real symmetric A. The resulting p' can then be written as the 
integral 

and the covariance matrix Y p of p then changes according to 
r p T p/ = T p + A 

This process is always allowed since adding a positive matrix A to a covariance 
matrix V gives a proper covariance matrix P\ For proofs of these statements, 
please see the proof of Theorem 6.1 and Lemma 6.2. 

Measurements The measurements on parts of multi-mode Gaussian states 
resulting again in a Gaussian state are exactly those, which can be described as 
projections on other Gaussian states. We will exploit this a bit when calculating 
the Schur complement in the next section. In the following we will discuss 
homodyne measurements. 
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4.2 Homodyne Measurements 

One of the experimentally well realisable measurements is the homodyne de- 
tection, where one mode of the measured N— mode state p is coupled to a local 
oscillator mode and measured together. The local oscillator is usually prepared 
in a coherent state with state vector \a p ) and the annihilation operator belong- 
ing to the local oscillator will be denoted by a p (with index p for pump mode). 




Figure 4.1: Homodyne measurement box 



With a homodyne measurement setup we are able to measure the expectation 
values of the quadratures of the measured mode 

<4>) = T2 • 

We present a small calculation without taking imperfect detectors into account. 
Behind the beam splitter we find the annihilation and creation operators of the 
modes b and b p which are composed of the incoming modes a belonging to p 
and dp belonging to the local oscillator with the transmission and reflection 
coefficients playing the role of the weighting prefactors. We get: 

bp = T* dp — R* a, 
b = Td + Rdp, 

with the commutator relations [b\b<] = |T| 2 + \R\ 2 when assuming [a, a p ] = 
[d,d P ] = [a) ,d p ] = [a) ,d P ] = 0. 




Figure 4.2: Beam splitter 

To set \T\ 2 + \R\ 2 = 1 therefore gives correct annihilation and creation opera- 
tors and is an energy conservation restriction for a non-absorbing mirror. The 
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expectation values of the number operators after passing the beam splitter are: 

(ft g ) = |T| 2 (ra a ) + \R\ 2 \a p \ 2 + \TRa p \ ■ (fie^ + aV^), 
(n B ) = |^| 2 (n a ) + |T| 2 |a P | 2 - \TRa p \ ■ (ae^ + aV^), 

with Q p the coherent state the pump mode a p was prepared in, and (j) a com- 
bination, <j) = (p T — (p R — <j) ap , of the phases of the transmission and reflection 
coefficient and the complex number a p respectively. If we first measure the 
number of photons coming out in the two modes behind the beam spitter and 
then subtract them from each other we find the number difference to be: 

(n~ b ) - (n s ) = 2\TRa p \(te i * + aV^) 

= V2\TRa p m<P)) = ^(x{®) : 

with, e.g., \T\ = \R\ = -±j. 

For a more accurate calculation involving imperfect photon detectors see Ref. 
[41], where it is shown that even in this case one can do precise measurements 
when using a strong local oscillator, e.g., a p is large. The output state and 
classical information after the homodyne detection is hence a projection of the 
iV-mode state on the localised states of the measured mode and the output is 
a Gaussian (N — l)-mode state. To go on we will need the covariance matrix 
and the displacement vectors of the, about a position x, localised state vectors. 



4.2.1 Moments of the Localised States 

We define the state vectors \ip x ,e) localised at x with width e. 



-{x-x'f 



<x'|^, e ) := : VJl / , (4.1) 



cxp 

V2vre 2 
with 

'I /'!„/. \|2 _ 1 



/r ' 2e0T 

So \ipxe) is n °t normalised but we find 



( dx'{x'\ip Xje ) = l. 
Jtr, 



Let us now calculate the first and second moments of the state vector \ipi x 

1 

tr[\lp Xj e){i> x ,e\] 



(%) = , ri ; TTi n / Ax' (x'\x\ip X)e )(^x,e\ x ') 



le^pK j dx'x 



R 

exp(^^ 



R 27T6 2 



— 1= f dx" (x" + x) exp ( X ) 
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because of the symmetry of the Gaussian. 



<x 2 > 



2e\pK / dx'x' 2 

J R 



cxp 



-(z-x') 2 



2vre 2 



1 



C0T 

1 



2 y — 5 — ) 



+ X 2 ey/7T 



2 e 



tr[\tJj Xie )(iJ} X£ \] 



dx {x'\p\lp X)e } {i>x,e\x') 



R 



2e^/vF y dx' \ -i 



d{x'\i) a 
dx 1 



{i>x,e\ x ') 



R 



dx' (x-x') \{i> x , e \x')\ 2 = 0, 



(p 



-:2 • 



2e^ / dx' ((-i) 



dx' 2 



2^ [ ^(^Ix-jfS ? 



R 



1 

2? ! 



= 2ie^ir dx — — ! (^eX) 

V dx' 



ie0r / dx'Ka/IV'x.e)! = « 



1 * 



¥1 



2e0r 2 



— ^ 

(px) = (xp)* = —. 



The displacement vector and the covariance matrix of the state vector 
are: 



(x) \ _ f x 




2((x 2 )-(x) 2 ) 2(xp)-i \ _ ( e 2 
2{px)+i 2((p 2 )- (p) 2 ) ) ~ { 1 
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This result is very convincing since for a well localised state vector \tp x e ) with 
e — > 0, the mean value should be x and the variance of the position should be 
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small. On the other hand we expect the variance of the momentum to grow, 
following the uncertainty relation. The uncertainty relation can be expressed 
in terms of the covariance matrix : det7 > 1 which is in our case fulfilled 

2 

for every e : detj Xi£ = % = 1. To go on we will use the covariance matrix 



t _ I ^ \ - r ( 



r i]£ = a 7 Xj£ ( q 2 ) an d the displacement D X)E = a d 

In the real world we can not do precise position measurements since our detec- 
tors have imperfect efficiency rates and always some dark counts. If we assume 
that in the limit of large numbers the measurement outcomes are approximately 
Gaussian distributed the following argument is true, even when the width e of 
the distribution is nonzero. 



Schur Complement 

We may now calculate what the outcoming state's covariance matrix is, after 
realising a projection on an arbitrary Gaussian one-mode state. For pab a 

A C 



Gaussian two-mode state with covariance matrix r = I ^ ) and dis 

placement D we have 



X P (U,Cb) := tr[W iPAB \ = exp + iD T tj , with £ T = e R 4 . 

Similarly we have for the Gaussian one-mode state u> on whom we project a 
covariance matrix T w and displacement D w 

Xu (vb) := tr[W v uj\ = exp ( -fl^l) , with r] T = V T B € H 4 .. 



With the expansion in Eq. (2.22) the state p' can be written as: 



P = tr B [p AB (t®u> B )} 

1 I' A t . (± 



tr B 
1 



J d 4 e xpi-Ofy (i® d 2 VB xU-vb)w vb \ 



far? 

l — f d^f d 2 VB Xp(-$)Xu,(-VB)Wt A tr B [Ws B W VB ] 



and with the trace formula for the Weyl operators in Lemma 2.3 we have 

p' = ^ / dU (J^ j d£ B x P (-0 xUte)) w u . 
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We find the characteristic function of the output state p' to be: 
Xp'(-&l) = J d£ B Xp(~0 Xw(£b) 

1 (a* ( ? V Z\ ( $?»Zb 
' d^2 exp — exp 1 



2?r ./ "V 4 / V 

x exp ^-il3 T £) exp (+iDl£B 

= 2^ exp 1 4 *-°a6i J / d £s exp ' 

xexp (-S^)exp {-i{pl-Dl)iB 

Here we realise that only projections on Gaussian states can transform all Gaus- 
sian states again to Gaussian states. If the characteristic function of to had third 
powers of the variable £ in its exponential, the characteristic function of p' could 
take a non-Gaussian form as well, for some Gaussian states p. We go on and 
transform the integration variables £b f] + C^ 1 ^ and include a transforma- 
tion constant J\f and get: 



Xp'{-U) = 7^ exp I — ^ zD A 6l I exp I I 

x exp exp f-i(D| - 



x y d^s exp 
x exp 



2 

'^cf- u, (r w + s)(7- 1 c fl 



4 

7? r (r, + b) c-^b 



exp (-^)exp {-i{D T B - Dl)C-\ B 



and set r/ = — [T^ + B] 1 \ C t ^a + 2i(D B — D w ) ) so that linear terms of ^ B in 
the integral vanish and the integral shortens to 



d£ B exp 



e B C- 1T (T w + B) C'^B 



with / being a real number. Using that B and T w are symmetric we finally 
get: 

Xp'(-U) = ^ ex P( 4 zL> a6i I 

x exp r ^Mk+^L^ ^ exp ^Tc[T u + B]- 1 ^ - 4 
x exp (- (d t b - DA \T U + B)-\D B - D u 
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The £^4— independent factors in xp' have to factor to one since the normalisation 
condition gives Xp'(O) = 1- Hence 



X p '(£a) = exp (- ^^ + iD^jy) exp 



x exp (-i&C[r w + B]-\D B - D u 
We find the displacement of the state after the projection to be 

D' = D A -C(T w + Br l (D B -D^) (4.2) 
and the covariance matrix 

r' = a - c (r w + By 1 c T . (4.3) 

The new V is a combination of submatrices of V and is known as the Schur 
complement of the matrix B + in f = T + ^ ~\ , see Ref . [6] . We 

will need it in the next chapter when discussing about measuring the degree of 
entanglement of an unknown state. 

Example 8 : Let us calculate V for the special covariance matrix 
r r e = ( ^> ) of the localised states. 



e 2 / 

We find with B = [ ^ ? 2 ) the inverse 



(T^ + B)- 1 



b 2 h 

bx + h b 2 



b 2 b 3 + e 2 
6 3£ 2 + e 4 _ he 2 

b 2 e 2 he 2 + 1 



(6 ie 2 + 1 )(5 3 + e 2)_ e 2 6 2- 

We finally get the covariance matrix after the projection 

(b 3 e 2 + e* -b 2 e 2 \ 
r , _ , r V he 2 + l) T 

(^e 2 + l)(b 3 + e 2 ) - e 2 b 2 ' 

Since the T x ^ e do not really depend on x, V is the same for every 
projection on any localised state, while the displacements are 
different for different projections. For e — > the above inverse 
becomes the Moore-Penrose inverse, generalising matrix inversion 
to singular matrices. 
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4.2.2 Projection on a Coherent State by Homodyne Measure- 
ment 

In the next chapter we will need an experimentally realisable setup to project 
one mode of a Gaussian multimode photon state on a coherent state. To do 
this projection we will use homodyne detection because experimentally it can 
often be achieved perfectly while other projection strategies are not as efficient. 
We discuss the case where the state p to be measured has two modes and 
we want to project it on a coherent state of its second mode. To realise this 
projection with homodyne measurements we have to find a scheme allowing us 
to do so. We will show that the following setup is convenient. 




Pout 



Figure 4.3: Measurement setup 

The mathematical description of this processing is done by unitaries followed 
by homodyne measurements in the second and third mode. 

&y,e\c ftMB Uc U§C (pAB ® |0)(0| C ) U^c tj}f \^ x>e ) B \^ e ) C . 

With all of these ingredients being Gaussian we can identify the outcoming state 
by using only the covariance matrix representation of all operations. We start 

A C 



with a two mode state p having an arbitrary covariance matrix T 



C T B 



and displacement D 




the vacuum 1 0> <0| with CM T = 1 and 



vanishing displacement (see Section 3.1). We will also need the localised states, 
with width e, realised by the homodyne measurements all having the covariance 

matrix T e = ^ <^ e 2 ^) an< ^ different displacements D x = f ^ ~\ . We apply 
a 50 : 50 beam splitter described by the symplectic transformation 



S BS 

a one mode ^— phase shifter 



J_ / 1 -1 
V2 V 1 1 



p 



si 

2 



1 
-1 



a 



and send two modes to homodyne measurement boxes. We calculate step by 
step what happens on the covariance matrix and the displacement of p passing 
this setup. 
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We abbreviate 
M 



M 1 



2 

1 + 6 3 + 2e 2 



62 1 + h + 2 



6 2 

1 + 7 



det 5 + 63(1 + 4 ) + + 2e 2 ) + 2(e + ±) 2 + 1 ' 



And get for 
and for I^ 4 ) we find: 

r( 4 ) = 



D A -V2CM -^D B -D X 



j=aD B - a{B - l)M 



T2 DB ~ D » 



A-CMC T c_{t-M{B- l))a q 

C T <r[{B+l)-(B-l)M{B-l)\a 
2 



As the last step we implement a homodyne measurement on the former third 
mode with the outcome y. The final displacement and covariance matrix are 
given by 



2}(4) ^ £(5) = g A _ CMy/2 



1 

V2 



C_ 

V2 



[1 — M(B — l))a 7 



a [(B + t)-(B- l)M(B - 1)] a 1 



x(-^aD B - cr(B - 1)M 
V2 



Dy), 



r W 1 — > = A-CMC T - -^{1- M{B -l))a T 

v2 

a [(B + 1) - (B - t)M(B - 1)] a 1 



xa(t-(B-t)M)^= 



Surprisingly this expression reduces to 

r( 5 ) = a-c[b + iy 1 c t , 



(4.4) 



where no e appear anymore, meaning that even a huge width of a Gaussian 
distributed measurement outcome is not relevant for the resulting covariance 
matrix. From the covariance matrix and the displacement of the resulting state 
p' we can read off that it would have been the same to project the initial state 
p on coherent state vectors \a(x,y,e)). The complex number a depends on 
the outcomes of the homodyne measurements: x and y, and on the Gaussian 
detection distribution with width e, with all those quantities appearing in 



Chapter 5 



Measuring Entanglement of 
Two-Mode States 



To determine whether a given unknown Gaussian two-mode state is entangled 
or not and how much, one can measure all entries of the covariance matrix, 
i.e., ten real numbers. But since we only want to know a special property 
of the state it could be that fewer measurements are necessary to answer the 
same question. We propose a scheme which makes it possible to determine the 
symplectic eigenvalues of the partially transposed covariance matrix with nine 
measurements. In the two-mode case, where the PPT criterion is necessary and 
sufficient, these symplectic eigenvalues are adequate to find out the degree of 
entanglement, measured in the logarithmic negativity. 

5.1 Measurement Scheme 

Given an unknown covariance matrix 1 7 of a Gaussian two-mode state 7 = 
A C \ 

C T B I w ^ ^' ^ rea ^ s y mme t r ic an d C real 2x2 matrices. 

Step 1 Measure all entries of A = ( 0,1 ° 2 ] an d determine a = V det A = 

V a 2 a 3 ) 

a\a-$ — a\. Experimentally this can be done by simple position and momentum 
measurements with homodyne detection. The measurement outcomes of a posi- 
tion/momentum measurement will approximately be Gaussian distributed. The 
variance of the measured distribution is an estimator for the covariance matrix 
element a\/a^. To measure the off-diagonal element 02 we apply a j— phase 

shifter = -4= ( , | so that the diagonal element An transforms to a 
4 V2 v l l J 

linear combination of all elements of A namely A' u = a\ + 03 — 2a<i- When 
measuring the position of the transformed state one will approximately get a 
Gaussian distribution with variance ^4' n from which we can determine the value 
of a 2 . 



We will use again the small covariance matrices 7. Keep in mind that the connection 



between the capital and the small CM is given by F = aja T . 
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Step 2 Measure all entries of B 
hb 3 - b\. 



, 1 2 ) and determine 6 = \/det B 
»2 b 3 



Step 3 Calculate and implement the local symplectic transformation which 
brings A and B to their diagonal form. 



S = Sa © Sb ■ 7 i— ► 7' = 5 , 7<S 



al C" 
C' T 61 



with C" = SaCS^. The experimental setup depends on the matrices A and B 
but the transformations Sa and Sb are generally phase shifters with variable 
cos 6 — sin t 



phase 



P m 



sm< 



COS i 



which should be available for every 



Step 4 As we saw in Subsection 4.2.2 it is possible to project a two (or more) 
mode system on coherent states of its second mode using homodyne measure- 
ments. Although the projections done with the proposed setup change with the 
measurement outcomes of the homodyne boxes, it is always a projection on a 
coherent state. The CM of the measured state changes in every case according 
to 

y ^ 7 " = al _( 6+ i)-i C 'c /T . 

giving a proper one-mode covariance matrix. Apply the projection physically. 

Remark: The different displacements of the coherent states have no effect on 
the covariance matrix but only on the displacement the resulting state has. 



Step 5 Measure all entries of 7" = ( ^ 1 ^ 2 ) . Since 7" is symmetric we only 

V 92 93 J 

have to measure three times, one less as when measuring the non-symmetric 
matrix C. Calculate the matrix C'C lT = (b + 1)(— 7" + al) and the absolute 
value of det C: 



1 /t/T\ 



det [C'C 



(detC) 2 = (det S A CSl) 2 

(det Cf = {b+ lf[{a -gi)(a- g 3 ) - g 2 ]. 



Step 6 With the following calculation we determine the determinant of the 
initial 7: 



det 



det 



1 

1 
al 



al a 
C' T 61 



1 



\aa T 

61 



6 2 det 



1 



al - f C'C lT 
6 



det [(6 + 1)7" - al] = det 7' = det 7. 
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Step 7 We now know the Simon normal form of the matrix 7 namely 



SNF 



( a 





c 


\ 





a 





d 


c 





& 





V 


d 





b ) 



1 



with a from Step 1, 6 (2) and c,d determined by the equations cd = detC (5) 
and det 7 = a 2 b 2 + c 2 d 2 — ab(c 2 + d 2 ) (6). There is a little ambiguity in 5 since 
the sign of the determinant of C is not fixed. 



5.2 Determining the Degree of Entanglement 

We remember that every two-mode covariance matrix can be brought to Simon 
normal form only by local symplectic transformations which do not change the 
degree of entanglement. The symplectic eigenvalues of the covariance matrix in 
Simon normal form are 



(7 )l,2 



-j= \Ja 2 + b 2 + 2cd ± yft^ 2 + b 2 + 2cdf 



4 det 7. 



(5.1) 



To decide whether the given state 7 is entangled or not we use the PPT-criterion 
stating that an entangled Gaussian state must have a partial transpose which 
violates the uncertainty relation for covariance matrices. The symplectic eigen- 
values of the partially transposed CM can be calculated using again the Simon 
invariants via : 
1 

V2 



h SNF )Ii 



— A/a 2 + b 2 - 2cd ± v / (a 2 + 6 2 -2cd) 2 -4det7. 



(5.2) 



From the preceding steps we have all ingredients to determine the values of 
(j snf )i j 2 an d ('~Y SNF )'i J 2 with the little subtlety that we do not know the sign 
of cd. We first take the positive sign for det C and calculate all eigenvalues. If 
all eigenvalues in Eq. (5.1) and (5.2) are greater than one, the state was not 
entangled. This does not change when taking the negative sign of det C since 
the eigenvalues of (7 SAr ^) TA and j SNF then just interchange. We hence know 
that the state is not entangled but can not decide if the state has Simon normal 
form 



SNFi 



fa 





c 







a 
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c 





b 
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d 





b J 



or 7 



SNF 2 
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a 





— c 


M 







a 





d 




— c 





b 





\ 





d 





b J 



If one of the eigenvalues was smaller than one then this eigenvalue belongs to 
(j snf } Ta and the state is entangled. In fact we are able to determine the 
degree of entanglement the state posses with only nine measurements instead 
of ten when measuring the whole state. The exact degree of entanglement can 
be calculated with Eq. (3.17) 

EM(l) =-^min(ln 7j TA ,0). 

i 

and the symplectic eigenvalues from above. 
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5.3 Discussion of the Measurement Strategies 

To experimentally realise a measurement is often a really expensive adventure, 
but the costs could be reduced when performing less measurements. The easiest 
way to determine the symplectic eigenvalues of a covariance matrix would be 
to measure just all elements of the CM; that makes ten different kinds of mea- 
surements. But why not learn from the available results of previous measure- 
ments? As we showed in the previous steps, it is possible to get the symplectic 
eigenvalues with less queries, when adjusting the strategy dependent on the in- 
formation we already gained. For future work, a challenging task would be the 
optimisation of such strategies. One could even try to proof how many kinds 
of measurements one necessarily has to perform to determine the symplectic 
eigenvalues of the covariance matrix of a given state, when it is allowed to learn 
during the process. 

But maybe our strategy is not that good because the variances of the measured 
covariance matrix entries could be worse as in the ten number case. It could be, 
that when measuring all entries of the CM one has to measure every entry Nio 
times to get the same results and the same variances as with our nine number 
idea for Ng >> N\q measurements. The energy saved when only measuring 
nine instead of ten entries is then spend on more tries for every element. 

We simulated both strategies, assuming that 7 is the matrix 



7 = 



/ 3.5 





2.5 










3 





-2.5 




2.5 





3.5 







V 


-2.5 





3 


/ 



The calculated symplectic eigenvalues of the given 7 are 

71 = 2.345 72 = 1.732 (5.3) 
and the symplectic eigenvalues of j Ta are given by 

7 p = 0.707 -Jz A = 5 - 745 > ( 5 - 4 ) 

thus 7 is entangled. To determine the entries of the covariance matrix we sim- 
ulated N = 10 2 , 10 3 , 10 4 , 10 5 , 10 6 measurements in total. For the usual ten 
number strategy we measured every entry N\q = yg times, while for the nine 
number strategy we measured every single quantity Ng = ^ times. The es- 
timations of the entries of the covariance matrix where done using ten (nine) 
Gaussian probability distributions, centered about zero. The variances of the 
first four Gaussians where the diagonal entries of the covariance matrix 7, which 
one can experimentally determine by position or momentum measurements. For 
the off-diagonal entries one first has to apply phase shifters, to bring them on 
the diagonal. Measuring again the position of such a transformed state gave 
values for those entries. For every expectation value we took the average over 
all N measurements and calculated the variance of the measurement outcomes. 
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These results are the estimators for the covariance matrix entries. The whole 
procedure was done M times for M = 1000 measurements. From the estimated 
covariance matrix we calculated the symplectic eigenvalues of 7 and r y TA , each 
M times, using Eq. (5.1) and Eq. (5.2) respectively and got an estimator for 
each of them. We compared the average of these measurement outcomes to the 
exact values and found the following result. 

With the usual ten measurements the estimation of the smallest symplectic 
eigenvalue Ji A , from which one can determine the degree of entanglement of 
the Gaussian state, was a bit better. The decrease of the abberation of the 
estimated symplectic eigenvalue in comparison to the exact value was approx- 
imate equal in both strategies. The standard deviation of the measured 7^ 
from the exact one, is calculated using the M estimated values of 7 1 A 




2 3 4 5 6 log 10 iV 



Figure 5.1: Ten Measurements versus Nine 

rp 

The figure shows the standard deviation A 7-^ for the ten measurement strategy 
(lower curve) and the nine measurement strategy (upper curve). The behaviour 
of the variances for the total number of measurements N = 10 2 , 10 3 , 10 , 10 5 , 10 6 
is shown on a logarithmic scale. Unfortunately our nine measurement strategy 
gives worse results than the ten measurement strategy. With increasing N the 
variances of the symplectic eigenvalue y^ A estimated with the two strategies 
decreases, as it should be. 
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5.4 Collective Measurements 

Finally we note, that it is by no means clear that single mode measurements, 
as we employed in the proposed setup, turn out to give better measurement 
strategies than a measurement using more copies of a state. Sometimes it can 
be useful to implement collective measurements, meaning that not only one 
copy of a state is measured but two or more copies are merged and processed 
together. The measurements of those multimode states could allow to further 
reduce the number of queries. We often want to know only a special quantity 
composed of entries of the covariance matrix. In this case, with collective 
measurements it is indeed possible to reduce the number of measurements, as 
the following example shows. 

Example 9 : Let 7 be a two-mode covariance matrix of the form 

1= {c T Bj- 

We are interested in the determinant of A + B and would naively 
measure the three entries of each matrix A and B. Employing collec- 
tive measurements instead could give the desired determinant with 
only three measurements, as the following strategy shows. Take two 
copies of a state with covariance matrix 7 and apply a beam splitter 
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We can just measure the three entries of the principal matrix B + A, 
where for the offdiagonal entry of B + A one has to implement an 
additional j— phase shifter on the second mode. With these three 
measurements we can calculate the determinant without knowing 
what the independent values of A and B are. We hence saved half 
of the costs with just a little trick. 



Chapter 6 
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We explain how separability can be formulated on covariance matrix level and 
introduce the advantageous concept of entanglement witnesses (EW). We pro- 
pose a scheme to efficiently estimate the degree of entanglement of an unknown 
state by realisable experimental measurements. Global reference is the script 
on Gaussian states, Ref. [42], to be published. 

6.1 Separability 

We review some basic properties of covariance matrices and formulate another 
separability criterion on the CMs. First, we recall the definition of separability 
of states and extend it to M parties. 

Definition 6.1 (Separability) A state p is separable with respect to M par- 
ties iff it can be written (or approximated, e.g., in trace norm), with probabilities 
Pi > and J2i Pi = 1 an d proper density matrices pf^ belonging to the jth party, 
as 



that is, the closed convex hull of the M -product states. 

Fortunately the separability of states can be formulated on the level of their 
covariance matrices with the following theorem, proved in Ref. [45]. 

Theorem 6.1 (Separability of CMs) Let 7 be the covariance matrix of a 
state p, which is separable with respect to M parties. Then there exist proper 
covariance matrices 7W , . . . , 7^) corresponding to the M parties, such that 



Conversely, if this condition is satisfied, then the Gaussian state with covariance 
matrix 7 is separable. If the stated relation is fulfilled by a CM we will name 
the covariance matrix itself separable with respect to the M parties. 




7 > 7 1 ' © ... © 7 



(M) 
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Proof : For the first statement let the state p be decomposed p = Y2i Pi P % where 
all p 1 are M-product states and with probabilities pi > 0, ^2iPi = 1. For 7 
covariance matrix and d displacement vector of the state p it is 

tr[p{R k - {R k ),Ri - (Ri)}+) = lki + 2d k d x = J^pMi + 244), 

i 

where the 7* are the block diagonal covariance matrices (Lemma 3.2) 
7* = Qjj-ii'j 1 )^ and d l the first moments of the p % . 

The difference A ki = j kl - YliPHkl = 2 J2iPi d \ d \ - 2d k d\ is a positive definite 
matrix, since for every v 6 H 2N it is 



(V, Av) = 2 ^ Pi ^2 V k d k d \ V l ~ 2 ^2 V ^ d kdlVl 
i kl kl 

= 2j2 P i\(vJ)\ 2 -2\(v,d)\ 2 

i 

and with d = ^2iPid l we get 



(v, Av) = 2j2pi\(vJ)\ 2 -2 

i 

> 2£p i |(t? >t ?)| 2 -2£p?|(t?,d 

i i 

= 2j>|(£, ( f)| 2 (l-;p i )>0 



Thus it is indeed possible to find a direct sum of proper covariance matrices 
W e . . . 7W, namely the 7 (i) := (Y^iPil*)^, tnat fuim the stated relation 



7 

since 
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Pi 
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In order to proof the second part of the theorem let 7 be a covariance matrix 

M 

7 = 7 (i) + A , 
i=i 

with an arbitrary positive matrix A and A a (classical) normalised, about zero 
centred Gaussian probability distribution in H 2N with covariance matrix — : 

let u> be a density operator describing a Gaussian state with block diagonal 
covariance matrix ©jli 7^ an d vanishing displacement. The expansion of to 
then reads 
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where we have used the relation T = a^a T . We construct the state p with A 
and lo : 

d 2 ^ \( V )W}ujW v , 

and find that it is a Gaussian state with covariance matrix 7 since its charac- 
teristic function is of the form : 

X P (0 = tr\ P W i ] = J d 2 ^ X^triW^W^} 
= tr[uW ( ] [ d 2N q \{r))e- iriT ^ 



= e 4 e 4 = e 4 

We now have to proof that is separable with respect to the M parties. We 
will need the following lemma to complete the proof: 

Lemma 6.1 Let C be a closed convex set, p be a C valued function and g be 
a continuous strictly positive function g : IR, — > R + with dtg(t) = 1. Then 
the integral 

r+oo 

dtg(t)p(t) (6.1) 
is again an element of C . 

Proof : Let G be the integral of g, that is G' = g, which is strictly monotonically 
increasing and G : R — > [0, 1]. We observe that the integral Eq. (6.1) can be 
substituted to the form 

+00 r+00 rl 

dtg{t)p{t)= dtG\t)p{t)= dxpiG-^x)). 
00 J ~oo Jo 

By definition of the Riemann integral that is equal to 

N 

lim Y /P (G- 1 (x l ))Ax u (6.2) 

with = x\ < X2 < ... xjy + i = 1 and Axi = Xi + \ — Xi > and hence 
Y2i=i Axi = 1. The sum in Eq. (6.2) is a convex combination and thus an 
element of C. Since C is closed the limit is reached in C. □ 

In our problem the multi-dimensional function A can be brought to product 
form by orthogonal diagonalization. Its factors are then indeed continuous, 
normalised and strictly positive functions. Furthermore all u> v := W\ uj W rj are 
product states, because the covariance matrix of uj is of block diagonal form, 
and hence separable, while the operators only displace u without affecting 
its covariance matrix. We use the previous lemma iteratively and find that the 
state p = J d 2I %] A (7/) lo v is separable. This completes the proof of the theorem. 
□ 
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Definition 6.2 We define the closed convex sets 

T(R 2Ar ) = Set of covariance matrices j, elements of Sym(2N,]R J ) with 
the properties 7 > and 7 ± io > 0. 

^a\b(^ 2N ) = °f separable covariance matrices 7 G r(R 2Ar ) ( with 

respect to the split A\B), with 7 > 7a © 7_B f or some 
j A € r(R 2Jv ^),7s G r(R 2A ^) with N = N A + N B . 

Definition 6.3 The principal submatrices A r of a square matrix A are the 
matrices which result when one neglects one or more columns and rows always 
of the same number. 

Lemma 6.2 Let 7 be a covariance matrix, element of T(IR 2Ar ) (r A | B (R 2JV )) , 
P a positive matrix and a > 1 then 

1. y = 7 + p 

2. 7 = 07 

are gfe'M elements ofT(R 2N ) (T A \ B (R 2N )) and 

3. all principal submatrices 0/7 containing N r modes 
are elements ofT(U 2N -) (F A \ B (K 2Nr )). 

Proof : 1. Is obvious. 

2. 7" is still a real positive matrix. The uncertainty relation is fulfilled for a > 1 
since a^f + ia = (a — 1)7 + (7 + icr) > 0. If 7 was separable with 7 > 7a © Jb, 
7" is still separable since 7" > a^y A © "7_B with proper covariance matrices a^ A 
and 075. 

3. All principal submatrices of a positive matrix A are always positive since for 
all v G ]R 2Ar : > 0. Especially for all v r in the reduced spaces R, r 

with iV r < AT. Then < (w r |^4|v r ) = (v r \A r \v r ) where the A r are the principal 
submatrices of A depending on the reduction one chose. That applies for 7,. 
and also for 7 r + ia r , so that j r is again a CM on r(]R, 2Arr ). Additionally the 
positive matrix 7 — 7^ © Jb stays positive for reductions on N r modes. □ 



6.2 Entanglement Witnesses 

A relatively new method to distinguish separable from entangled states results 
from the insight that the set of covariance matrices and the set of separable CMs 
form closed convex subsets of the space of real symmetric 2A^ x 2A^ matrices. 
The convexity allows to describe these subsets completely by the intersection of 
closed half-spaces or equivalently by a family of linear inequalities. One family 
of such inequalities is given by the so called entanglement witnesses. They were 
first proposed in Ref. [43], further investigations can be found in Ref. [44] on 
discrete states and on continuous states in Ref. [42] . We consider bipartite CVS 
systems and start with the definition of appropriate entanglement witnesses. 
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Definition 6.4 We define the following sets of real symmetric 2N x2N -matrices 
Z(K 2N ) = {Z\ V 7 e T(R 27V ) : tr[jZ] > 1} 
2U| B (R 2JV ) = {Z|V 7 G r A | B (R 2JV ) : tr[ 7 Z] > 1,} 

where A\B denotes the split of a bipartite system with N = N A +N B modes and 
7 covariance matrix. All matrices Z £ ^^(IR 2 ^) will be called entanglement 
witnesses for the split A\B. 

The denned sets are again closed convex subsets with Z{Wi 2N ) C Z A]B (Il 2N ). 
Furthermore it is possible to assign to every entanglement witness a quadratic 
Hamiltonian. 

Definition 6.5 Let Z be an element of Z A \ B (M 2N ) and R = (x±, . . . , pn)- T 
We define the Hamiltonian 

N 

Z[R] : = 2 ZhlRkRl- 

k,l=l 

With the help of the trace we see how the matrix Z is related to the expectation 
value of the operator Z. Let 7 be the covariance matrix and d the displacement 
of a state p. The expectation value of any Z can be expressed in terms of the 
moments of the state: 

N 

{Z[R]) p = tr[Z[R]p] = 2 Z kl tr[R k R lP ] 

k,l=l 

N 

= Y Z kiHl + 2d k di) = tr[Zj] + 20? Zd. (6.3) 

k,l=l 

In the following lemma we see how entanglement witnesses detect entangled 
states. 

Lemma 6.3 For 7 covariance matrix of a bipartite system of N A + Ng = N 
modes, it is: 

7 is separable VZ G ^^(IR, 2 ^) : tr[Zj] > 1 

with respect to the split A\B. 

For Gaussian states only the covariance matrix determines the entanglement 
properties, thus a separable covariance matrix is necessary and sufficient for 
a Gaussian state p to be separable. For non-Gaussian states it is necessary 
to have a separable covariance matrix but not sufficient since higher moments 
could show entanglement. The lemma for non-Gaussian states is thus 

3Z e Z A \ B (K 2N ) : tr[Zj] < 1 7 is entangled 

=^ all CV states having the covariance 
matrix 7 are entangled. 

Proof : is clear from the definition of entanglement witnesses in Def. 6.4. 
In order to prove the <£= direction we will need an important theorem from the 
theory of convex sets which can be found in [46] given here without proof: 
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Theorem 6.2 (Separation of point from convex set) 

Let Y be a closed, convex, nonempty set of matrices, Y C Sy?n(2N, R) , let 
x G Sym(2N, B)\Y. Then there exists a matrix M G Sym(2N,H), M / 

tr[Afx] < min ir[Afy]. 

yeY 



The set of separable covariance matrices T^ B (JR, ) is such a closed convex set. 
In order to prove Lemma 6.3 we assume to have an entangled x and show that 
there has to exist aZ£ Z^ B (H 2N ) which then detects it and gives tr[Z x] < 1. 




Figure 6.1: Separation of point from convex set 
With Theorem 6.2 we find a matrix M G Sym(2N, R) so that 



tr\Mx] <c= min tr[M 7'] < tr[M j] V7 G r j4 , i3 (R : 

7'er A | S (R2iV) 



2A^ 



(6.4) 



Let us assume that the real number c is strictly positive. Then we get 
ir[M x] < 1 < tr[^ 7] V7r A | B (R 2Ar ) but if this is true M = Z has to be an 
element of -Z^|^(R 2Ar ), since it is a real symmetric matrix with the defining 
condition in Def. 6.4. We have found that if we take an x F A \ B (1R 2N ) 
we will find an Z G Z^| B (R 2Ar ) which gives tr[Z x] < 1 or equivalently 
G Z A \ B (K? N ) : tr[Zi\ > 1 => 7 G T j4 | i? (R 2 ^). 

We still have to discuss the case c < 0. Let us show that M is positive by 
contradiction. Suppose M were not positive and had an eigenvalue tjiq < 
with the corresponding eigenvector \i G R 2Ar . We construct the matrix 7I, = 
7fc^+A/ifc//j for an arbitrary A > and a separable matrix 7, which by Lemma 6.2 
is itself an element of T A | jB (R 2Ar ). We find that tr[M-y'] = tr[Mj] + Am /i 2 -> 
—00 for A — ► 00. With this construction we see that M cannot have any 
eigenvalue smaller than zero because Eq. (6.4) has to hold and tr[M *y'] = 
tr[Mj] + Am-o/i 2 would reach values smaller than tr[Mx] for sufficiently large 
A. Hence M is positive. 



lr The theorem holds for more general sets, we state only what we will need. 
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But then the following equation holds for every 7 G T^ B (H 2N ): 

tr[Mj] = tr[MS 1 WNF S T ]=tr[S T MS 1 WNF ] 
= tr[S T M S{ 1 WNF -l)]+tr[S T M S] 

> 0, 

since the first addend is again the trace of a (semi-)positive matrix, see Eq. 
(3.11) and Lemma 2.1, while the second addend is the trace of a positive matrix 
which is only zero if all eigenvalues of it were zero that is the matrix S T M S 
and hence the matrix M itself was the zero matrix. But Theorem 6.2 excluded 
the M = case. We conclude that there exists at least one nonzero M G 
Sym(2N, ft) so that for all 7 G T^ B (]R 2N ) we find a strictly positive c. This 
completes the proof of Lemma 6.3. □ 

Definition 6.6 The symplectic trace (denoted as str [A] ) of an even dimen- 
sional, positive, symmetric matrix A is defined as the sum of the symplectic 
eigenvalues (counted once). 

( A\ A2 \ 

Lemma 6.4 Let A = [ . T . I denote a real positive 2N x 2iV matrix, 

\ A 2 A 3 ) 

then the symplectic trace has the following properties: 

str[SAS T ] = str[A] VS G Sp{2N, R) invariance 
str[pA] = pstr[A] p > multiplicity 
str[A] < str[Ai] + str[A 3 ] 

Proof : The invariance is clear since symplectic transformations preserve all 
symplectic eigenvalues, additionally the symplectic eigenvalues just all multiply 
with the factor p so that multiplicity holds. The last property can be derived 
as a consequence of Theorem 6.3. □ 

But given a real symmetric 2iV x 2N matrix, is it an entanglement witness or 
not? The following theorem answers the question providing a check, based on 
the symplectic eigenvalues the given matrix possesses. 

Theorem 6.3 Let Z be a real symmetric 2N x 2iV matrix on a phase space of 
N = Na + Nb modes. Let Za and Zb denote the principal submatrices of Z 
belonging to the subsystems A and B respectively. Then : 

1. Z G Z(K 2N ) Z > and str[Z] > \ 



2. Z G Z A \b{^ 2N ) Z > and str[Z A ] + str[Z B ] > 1 



Note that str[Z] > | is stronger than str[ZA\ + str[B] > |. 
Proof : (following reference [42] ) 

1. =>■ : With the same argument as in the proof of Lemma 6.3 we find that Z G 
Z(Wl 2N ) has to be positive. Since Z > we can find a symplectic transformation 
to bring Z to Williamson normal form: 



Z WNF = SzZS T = Zo(B Q l {zi 0\ 



62 



CHAPTER 6. ENTANGLEMENT WITNESSES 



where Zq is a matrix with determinant zero and zero symplectic eigenvalues. We 
assume that the particular Z has a vanishing Zq part an go on. The symplectic 
trace of Z is now easily calculated: str[Z] = str[Z WNF ] = YliZi- Furthermore 
we have with the basis transformation R' i— ► R = S Z R! 

tr[Z[R'}p'} = tr[Z WNF [R]p] = 2Y,{{Ql) P + {Pf) P )z l 

i 

= 2^2(2(hi) p + l)zi > 2 str [Z], (6.5) 

i 

where hi = ajoj. The right hand side is the expectation value of the Hamiltonian 
of harmonic oscillators which is bounded from below by the ground state energy 
reached when p in Eq. (6.5) was the ground state. 

On the other hand the expectation value of Z is related to the state's moments 
via tr[Z[R']p'] = tr[Z-y'] + 2d! T Zd' > tr[Z-y'] > 1 for all 7' G r(IR 2Ar ), p' G 
S{TL) by Def. 6.4. Since the boundary in Eq. (6.5) is reached, we conclude the 
stated property of Z, namely 2str[Z] > 1. Also for all Z' identical to Z but 
with a nonzero Zq part it follows str[Z'] = str[Z] > |. 

1. <=: For all 7 G r(]R 2Ar ) and with Z WNF again Williamson normal form of 
Z, including a Zq part, we calculate: 

tr[Z 7 ] = tr[S T Z WNF S 1 ]=tr[ZQ 1 l Q] + Y,tr[zM\ 

i 

> ^^(Aj + Af), 

i 

where 7q and 7- are the principal submatrices of 7'. 70 has the size of Zq 
while all 7^ are 2x2 covariance matrices with the positive eigenvalues Xj and 
A 2 which fulfil Xj • A 2 > 1 for every submatrix since the determinant of 7^ is 
necessarily greater or equal to one. And since for every A > the relation 
< (y/X - ^) 2 = A + i - 2 holds, it follows that 

V 7 er(R 27V ) : tr[Zj] > ^ Zi (X\ + ^) > 2^ Zi = 2 str[Z] > 1. 

i 1 i 

Thus Z is an element of i?(]R 2Ar ) as stated. 

2. =^ : For Z G iJ^^lR 2 ^) we run the above argument (proof of Lemma 6.3) 
for the positivity. For a separable 7 we have the useful condition proved in 
Theorem 6.1: 7 > 7a © 7 s for some covariance matrices 7a and 7#- Thus 

tr[Zj} > tr[Z AlA ]+tr[Z BlB \ > 

and with Eq. (6.3) and pa, Pb the Gaussian states with covariance matrix 
1A-,1b respectively and vanishing displacement, follows 

tr[Zi\ = tr[Z A [R]p A ]+tr[Z B [R]p B }. 

Additionally it is Za,Z b > and hence irfZ/if^Jp^] > 2str[ZA\ and 
tr[Z B [R]p B ] > 2str[Z B ] from above, so we get 

tr [Zj\ > 2str[Z A ] + 2str[Z B ] 
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reached again when both subsystems are in ground state. By the precondition 
we also have tr\Z^\ > 1 for all 7 G ^a\b(^ 2N ) so that the condition on the 
symplectic trace of the entanglement witness is 

2str[Z A ] +2str[Z B ] > 1. 

2. 4= : We use the same idea as for 1. For all 7 G r j4 | B (R 27V ) and Zjf ^ 
and Z B VNF the symplectic diagonalized submatrices of Z it is : 

*r[Z 7 ] > tr[Z A1A }+tr[Z BlB ]=tr[Z A VNF 1 ' A )+tr[Z B VNF 1 ' B ] 



= tr[Z^ A ] + Y,tr[zttiC\ 

+tr[Z B 7 ' B ]+zZtr[zFH B ] 

i&B 

> 2j2 z f + 2 J2 z ? = 2 str \ z A] + 2 str[Z B ] > 1, 



finally completing the proof. 



□ 



We have now a powerful tool to decide whether a given matrix is an entangle- 
ment witness for the given split A\B or not. And with them we can restrict 
the sets a given covariance matrix is an element of. The picture shows how two 
entanglement witnesses Z\ and Z2 divide the set of covariance matrices in ]R, 2Ar 
into four subsets. If the trace of the given CM 7 together with the entanglement 



separable covariance matrices 
trace of Z±j is smaller than 1 
^ trace of Z27 is smaller than 1 




Figure 6.2: Entanglement witnesses dividing T(]R 2JV ) 



witness Z\ (Z2) was smaller than one, then 7 belongs to the red-lined (blue- 
lined) set respectively. If both measurements gave values smaller than one, it 
is in the squared set. In those cases all states having the covariance matrix 7 
are entangled. If, on the other hand, both entanglement witnesses gave values 
greater than one, we can not say whether the given matrix is entangled or not. 
But from the geometric representation we guess that optimal entanglement 
witnesses going exactly through the edge between separable with entangled 
CMs, exist and indeed they do. We will not introduce them here but only give 
references [47]. 



64 



CHAPTER 6. ENTANGLEMENT WITNESSES 



6.3 Estimation of the Degree of Entanglement 

We now come to the main results of this thesis. We have so far learned a 
lot about Gaussian states, their covariance matrices and their entanglement 
properties. When producing a Gaussian two-mode state in the lab usually one 
wants to make sure that the outcoming state behind the particular setup is 
really the desired state. Especially it is often necessary to check whether the 
state has the expected entanglement properties used for quantum cryptography 
or other purposes. For that reason one is interested in a simple measurement 
scheme telling the experimenter immediately what degree of entanglement the 
produced state has. In the third chapter we learned that for Gaussian states 
entanglement, quantified in the logarithmic negativity, can easily be calculated 
using only the covariance matrices of the states. So it should be possible to 
measure some properties of a CM and give an estimate of the degree of en- 
tanglement for the Gaussian states having this particular covariance matrix. 
In the previous section we introduced the concept of entanglement witnesses 
and saw that they are very practical and provide a natural way to distinguish 
separable from entangled states. But one could do even better. It is not only 
possible to decide whether a state is entangled or not but with practicable mea- 
surements one can estimate a lower bound of entanglement of a given unknown 
Gaussian state. We show how this can be done and start off with the definition 
of p— separability. 

Definition 6.7 For any 7 £ r(R 2Ar ) covariance matrix of a bipartite system 
A\B and any p E (0, 1] we define 



7 is p — separable 



— is separable 
p 



with respect to the split A\B respectively. The sets of p— separable covariance 
matrices are again closed, convex subsets of r(]R, 2Ar ) and will be denoted by 



The definition is sensible since 
covariance matrix properties, 
to I — separability and for p\ 



- with Lemma 6.2 and a = - 

p p 

Note, that by definition separability is equal 
> P2 follows: 7 is pi — separable 



> 1 still has 



7 IS P2 



separable and the sets T P MB (JS? N ) include each other like the shells of an onion 



r 1 

1 A\B 



(R : 



L 1 A\B 



(R : 



A\B 
2N\ -nO.2 
L 1 A\B 



(R 2N ) C r(K, 2JV ), as shown in the figure. 




(1-) separable CM 
0.5-separable CM 
O 0.2-separable CM 



Figure 6.3: p— Separability 
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Using the definition of p— separability Theorem 6.1 can be reformulated to 

Theorem 6.4 (p— Separability of CMs) Let 7 be a, with respect to the two 
parties A and B, p— separable covariance matrix, i.e. 7 6 r^ B (]R 2Ar ). Then 

there exist proper covariance matrices j A £ r(]R 2JV>1 ) and j B £ r(R 2Ars ) cor- 
responding to the parties, such that 

1>P {l A ®l A ). 

And finally, p— separability of covariance matrices can easily be expressed with 
entanglement witnesses as the following lemma shows: 

Lemma 6.5 For 7 CM of a bipartite system, p £ (0, 1] and Z being an entan- 
glement witnesses, we find the relation: 

7 isp- separable (A\B) VZ £ Z A \ B (B? N ) : tr[Zj] > p 

or 7 is not p - separable (A\B) -t=^ 3Z £ Z A \ B (K 2N ) : tr[Zj] < p. 

Proof : With the definition of p— separability, Theorem 6.1 and Lemma 6.3 the 
stated relations follow immediately. □ 

For p < 1, p— separability is weaker than the usual separability, but when blown 
up with a factor - they would reach the set of separable covariance matrices 
T^ B (M 2N ). In this sense, the smaller the number p is, the more entangled can 
the covariance matrix be. We observe, that when measuring an entanglement 
witness Z, the outcome p = tr\Z^j\ does not only distinguish separable (p > 1 
for all witnesses) from entangled (p < 1 for one particular witness) covariance 
matrices, furthermore it orders the entangled CMs in sets of p— separable CMs. 

An experimenter who produces Gaussian states with a certain degree of entan- 
glement is naturally curious if he succeeded in doing so. When he measures 
an entanglement witness Zq on a given unknown Gaussian state p having a 
vanishing displacement and covariance matrix 7 the outcome is given by 

m = (Z [R]) p = tr[Z j] <p:=m + e, Ve > 0. 

The measurement outcome tells him now immediately that the state is not 
p— separable for all e > 0. But maybe this is not the best estimate he could 
do. It is possible that with a different entanglement witness a smaller outcome 
m could be measured and the state would be even more entangled. Since he 
produced the Gaussian state p, it is not completely unknown to him; he actually 
has a strong conjecture what the state will be and can use his knowledge to 
find an minimal entanglement witness for the covariance matrix associated to 
the state p. In contrast to optimal entanglement witnesses it is not the witness 
going exactly through the edge between entangled and separable states, but 
it is the one that gives the best estimate of the degree of entanglement the 
supposed state has. Obviously the minimal entanglement witnesses depend on 
the covariance matrix for whom they are minimal. We will understand the 
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minimal entanglement witnesses Z m i n {pf) to a given bipartite CM 7 as being 
those witnesses Z m i n {^) G ^| B (lR, 2Ar ), which give 

tr[Z min (j)i] =p m in(l), 

while all other EWs Z £ Z A \ B (1R 2N ) give tr[Zj] > p m in- The number Pminil) 
is then the smallest value of the outcomes p, one could measure for the CM 7. 




Zopt 



Figure 6.4: Optimal and minimal entanglement witnesses 

In Fig. 6.4 the coloured areas are those where entanglement witnesses gave 
values smaller than one, so the CMs in these areas must be entangled. The 
minimal EW Z m i n {^ ) gave tr[Z min (j ) 70] = Pminilo) for the covariance matrix 
70, symbolised by the red dot. 

We see that it should be possible to give lower bounds for the degree of en- 
tanglement of a Gaussian state using entanglement witnesses. To estimate the 
degree of entanglement quantified in the logarithmic negativity, we will need 
that the PPT-criterion is necessary and sufficient, which is the case in 1 A x Nb 
mode systems. We learned in Lemma 3.6 that the logarithmic negativity of a 
Gaussian state can be calculated from its covariance matrix: 

EM{l) = - J^min(ln 7 ?\0). (6.6) 

i 

In the following lemma all the things we learned are merged together and we 
find how the outcomes of the measurements of entanglement witnesses bind the 
logarithmic negativity of the measured state. 

Lemma 6.6 In case the PPT-criterion is necessary and sufficient: 
If an experimental setup measuring the trace of an entanglement witness Z with 
a covariance matrix 7 belonging to an unknown Gaussian N— mode state gave 
the outcome m G (0, 1), the logarithmic negativity of the state is bounded from 
below by 

Exirt) > In-. 

m 

In the two-mode case, where only one symplectic eigenvalue of j Ta can be 
smaller than one, the minimal witness Z m i n giving the smallest possible value 
mmin gives exactly the logarithmic negativity the state possesses. 
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Proof : With the symplectic eigenvalues of the partially transposed covariance 
matrix j Ta defined in Eq. (3.15) and Lemma 6.5 we see that forall p £ (0, 1] 

Vi : 7 4 Ta > 1 <^=^ 7 is separable (PPT!), 

Vi : 7 4 Ta > p <J=^ 7 is p — separable •<=>• VZ G ^| B (lR 2Ar ) : tr[Zj] > p, 
3i : jJ A < p <J=^ 7 is not p — separable 3Z E Z^ B (]R 2N ) : tr[Zj] < p. 

From the measurement we know that there exists an entanglement witness 
giving the result m. We hence know that there exists at least one symplectic 
eigenvalue j[ A of j Ta which is smaller than p = m + e for all e > 0. Thus we 
can give a nontrivial lower bound for the logarithmic negativity of the state. 

JV 

Erfd) = ~ min(ln ^J A , 0) — min(ln j[ A , 0) 

i=2 

> — lnp = In 

m + e 

and since this inequality is true for p = m + e, Ve > 0, the proof of the first 
statement is complete. For the second part we assume a symplectic eigenvalue 
7^ A smaller than one and we get 

Ejvil) = — min(ln ^ A , 0) — min(ln 7^ , 0) 
= + ln^. 

Since the outcome of our measurement gave the value m m i n for the given co- 
variance matrix, we conclude that the symplectic eigenvalue is smaller than 
or equal to m m j n . On the other hand by definition of the minimal measurement 
outcome m min it is for all Z £ ^^(IR 2 ^) : tr[Zi\ > m min and hence for 
all symplectic eigenvalues r y i A > m m i n . Hence 7 : A is equal to m m in and the 
logarithmic negativity is given as stated. □ 

The lemma provides an easy method to estimate a lower bound of entangle- 
ment by measuring different entanglement witnesses. The estimated logarithmic 
negativity becomes larger the smaller the minimal measurement outcome was, 
respecting the ordering induced by the sets of p— separable covariance matri- 
ces. Entanglement witnesses also detect bound entangled states, being those 
entangled states having a positive partial transpose. If existent, they destroy 
the sufficiency of the PPT-criterion and measured in terms of the logarithmic 
negativity they have a vanishing degree of entanglement. Unless bound en- 
tangled states exist the logarithmic negativity is strictly positive for entangled 
covariance matrices. A possible goal for future research could be a computable 
entanglement measure also respecting the p— ordering but giving nonzero values 
for bound entangled states. 

Summarising our achievements, we see that an experimenter, wanting to know 
the degree of entanglement his produced Gaussian state possesses, simply has 
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to measure some entanglement witnesses and the results instantly tell him how 
much the state is minimally entangled. With little additional knowledge about 
the state it is furthermore possible to optimise the measurements, trying only 
minimal entanglement witnesses. In the following example we see that some- 
times even the knowledge of the structure of the covariance matrix is sufficient 
to determine minimal entanglement witnesses for it. 

Example 10 : Let 7 be of the special form 



7 



/ a 





b 





\ 





a 





-b 




b 





a 







\ 


-b 





a 





with a > 1 and < b < a. 



Then the symplectic eigenvalues of the ^ Ta are 7J" 4 = a — b and 

rp 

7 2 A = a + b. The logarithmic negativity of 7 is then given by 
Ej\f(j) = In ^£ and the minimal measurement outcome for an 
entanglement witness is m m i n = a — b. An entanglement witness 
which allows to measure the smallest value m m i n is given by: 





( 


1 





-1 


°\ 


1 







1 





1 


4 




-1 





1 







\ 





1 





1 / 



with tr\Z min ^\ 



A small check that Z m i n is really an entanglement witness in 



J A\B 



(ft 



27V \ 



Z u 



eigenvalue of Za = Zb 



is positive since its eigenvalues are with Eq. 

■1 

is z\ 



(3.7) given by Ai j3 = > and A 2 , 4 = M > 0. The symplectic 



-7 hence the condi- 



1 

1 ; is zi = * 

tion on the symplectic trace of the submatrices is fulfilled since 



str[ZA\ + str[Z B ] 

by 



1 



The Hamiltonian to be measured is given 



Z[R] 



2 ZkiRk Ri 



(x\ - x 2 ) 2 , (pi +h) 



+ 



fei=i 



and is just a beam splitter mixing the two modes followed by a 
measurement of the position and momentum. 



In the example we were sure to have a CM of the given form with entries a 
and b unknown to us. We only used the knowledge about the structure of the 
CM and found a possible minimal entanglement witness. When measuring the 
proposed EW we got the exact logarithmic negativity independent of the actual 
values in 7. Similar optimisations could be done for every particular structure 
a CM can have and the better the knowledge of the CM the more exact one 
can determine the appropriate minimal entanglement witness. 
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6.4 Duan Criterion 

The Duan-criterion proposed in 2000 by L.M. Duan et. al., Ref. [48], is a 
special case of the theory of entanglement witnesses. 

Lemma 6.7 (Duan-criterion) 

The criterion states that given a separable two-mode state p = ^2; pi pf (8) pf 
every pair of operators 

1 

u = \a\xi H — X2, 

° oeE,a/0, (6.7) 

v = \a\pi p 2 

a 

has to fulfil the relation 

A 2 u + A 2 v>a 2 + ij, (6.8) 
a z 

where A 2 u = (u 2 ) — (u) 2 and similarly A 2 v are evaluated in the state p. 
A necessary condition on a state to be entangled is thus that there exists a value 
a such, that the above relation is violated. For Gaussian states, Duan et. al. 
showed that the condition is also sufficient. 



(6.9) 

a 2 (xj) + -^(x 2 2 ) + a 2 {p 2 ) + + 2 ^ {{ Xl x 2 ) - (pip 2 )) ■ 



We translate Duan's proposal in the language of entanglement witnesses. With- 
out loss of generality we assume that the states to be measured have vanishing 
displacements. Then the measured quantity following Duan's proposal is given 
by 

A 2 u + A 2 v = (u 2 ) + (v 2 ) 

^{x 2 )+a 2 {p 2 ) + \{p 2 ) + 2^ 
a z a z a 

In the language of covariance matrices we find 

A 2 u + A 2 v = ^— (711 + 722) + 77-9 (733 + 744) + — (713 - 724) • (6.10) 
z 2a z a 

We put Eq. (6.9) in an appropriate form to read off the corresponding entan- 
glement witnesses. 

A 2 u + A 2 v 1 [ 2 , ,1 

n , 2 , n a v7n + 722) + — (733 + 744) 
Aa + . 

1 \a\ . 

2 — (713 - 724) 



1 < 



a? + 4 

a- 



+ 



2(a 2 + \ 



(6.11) 



One family of entanglement witnesses Z a which give tr[Z a p] 



A 2 u+A 2 v 

.2 1 1 
a 1 — 7 



when 



measured in an about zero centered state p has the matrix representation 







a 2 





H 

a 





\ 


1 







a 2 





_H 

a 




2(a 2 + 




M 

a 





1 









I 





_H 

a 
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/ 
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By construction these witnesses, with arbitrary a / 0, give for all separable 
states values greater or equal to one, see Eq. (6.8). For all a ^ Z a is then 
indeed an entanglement witness since it is a positive, symmetric matrix fulfilling 
the condition on the symplectic trace of its submatrices 

1 11 

str[Z aA ] + str[Z aB ] = (a 2 + — ) = - 

2(a 2 + jz) a 2 2 

and give values greater or equal to one for all separable covariance matrices. 
The sufficiency is a bit surprising since the family of these entanglement wit- 
nesses is not the whole set Za\r(JS? n ), for a proof see the original paper [48]. 

We conclude that the theory of entanglement witnesses is more general than 
the Duan-criterion, since it can be used for multi-mode states as well. We also 
note that the set of necessary EWs could be smaller than defined in Def. 6.4 
and characterised by Theorem 6.3, as the Duan-criterion shows. Further in- 
vestigations should analyse which entanglement witnesses are redundant. For 
example the set of EWs could be equipped with an equivalence relation collect- 
ing entanglement witnesses, which can be transformed into each other by local 
symplectic transformations in equivalence classes. This would not change the 
theory of EWs since the entanglement properties of the states are preserved 
under LOCC. 



Chapter 7 

Summary and Outlook 



In the thesis we investigated entanglement properties of Gaussian states and 
how well one can estimate them for a given state with real experimentally avail- 
able measurements. In the second chapter we started to dive into the language 
of continuous variable states. We learned what symplectic transformations are 
and how they are related to unitary transformations. We introduced the quan- 
tum characteristic function, in analogy to the classical Fourier transformed of 
a probability distribution, describing a quantum statistical system. Finally we 
formulated the definition of entanglement of bipartite states and cited the fa- 
mous PPT-criterion allowing to judge if a state is entangled or not. In the third 
chapter we defined Gaussian states and found that well known states like the 
coherent, the squeezed and the thermal states of a harmonic oscillator belong 
to the family of Gaussian states. We recognised that for Gaussian states the co- 
variance matrix plays an extraordinary role, describing all important properties 
like the Heisenberg uncertainty relation, purity, squeezing and especially entan- 
glement of the state. We investigated two normal forms of covariance matrices 
and formulated conditions on the CM's characteristic numbers connected to 
them. We went on to exploit general operations one could implement on Gaus- 
sian states, without destroying this important property. We saw that in fact 
many experimentally available operations are of Gaussian nature, and learned 
what homodyne measurements do and how a projection of a two mode Gaus- 
sian state on coherent states of its second mode changes the state's covariance 
matrix. 

In the fifth chapter we gave a first try on estimating the entanglement properties 
of an unknown Gaussian two-mode state. Instead of measuring all entries of the 
state's covariance matrix we proposed a scheme allowing to calculate the exact 
symplectic eigenvalues of the CM and its partial transposed with nine kinds 
of measurements. With them we can determine the degree of entanglement 
the Gaussian state possesses, quantified in the logarithmic negativity. Many 
possible ways lead from this ansatz; one could ask if it is possible to determine 
or well estimate the entanglement properties of a given Gaussian two-mode 
state with even less than nine measurements. Other properties like purity 
or squeezing could be investigated with similar strategies and probably less 
measurements. In the next step one should try to find schemes for entanglement 
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estimation for more mode states and non-Gaussians. 

Finally, we introduced the fantastic concept of entanglement witnesses allowing 
us to give a lower bound for the logarithmic negativity of a Gaussian state 
by measuring quadratic Hamiltonians. Entanglement witnesses show a lot of 
nice properties, e.g., they detect bound states, being those entangled states 
having a PPT. In comparison to other measurement requirements they are 
easily available and with a minimal entanglement witness for a given CM, it 
is actually possible to give the exact logarithmic negativity. Many questions 
concerning EWs remain open and should be answered in the near future. We 
give now some proposals for further research issues: 

When discussing the Duan-criterion in Section 6.4 we saw that the set of 
entanglement witnesses is bigger than necessary and should be reduced to 
a minimal set. A first step in this direction would be a equivalence relation 
grouping all EWs equivalent up to local symplectic transformations in 
equivalence classes. Further reductions can probably be found. 

Bound entangled states could maybe be characterised by the entanglement 
witnesses detecting them. Is it then possible to find a nontrivial lower and 
upper bound for the minimal Pmin 

for all bound entangled states? 

In our discussion we only used the logarithmic negativity. One could check 
if the entanglement witnesses allow statements about the degree of entan- 
glement of Gaussian states measured in other entanglement measures. 
Additionally the p— separability seems to imply an ordering concerning 
the degree of entanglement in the set of covariance matrices. It should 
be investigated if the minimal p m in{l) itself is a proper entanglement 
measure for the CM 7 or if known entanglement measures respect the 
p— ordering as the logarithmic negativity does. 

We gave a proof that for two-mode Gaussian states the logarithmic nega- 
tivity is exactly given by the logarithm of the inverse minimal p m in (7) one 
could measure for the CM 7 belonging to the state. Is that also possible 
for more mode states? And are there other strategies to determine the 
minimal entanglement witness for a particular class of covariance matrices 
like in Example 9? 

The entanglement witnesses were formulated only on the covariance ma- 
trices. Hence they decide about the entanglement properties of Gaussian 
states and give a sufficient condition for non-Gaussian states having an 
entangled CM. Is it possible to generalise the concept of entanglement wit- 
nesses to higher moments so that it applies to other classes of continuous 
variable states? 

And finally, for systems split in more than two parties what definition of 
entanglement shall one choose and how can entanglement witnesses used 
for multi-partite entanglement? 



Appendix A 

Definitions and Proofs 



Definition A.l The Poisson bracket is defined by 
' r-** dxi dpi dxi dpi 

i=l 

The commutator of two M x M matrices or operators A and B is defined by 
[A, B] = AB - BA. 

The anti- commutator of two M x M matrices or operators A and B is defined 
by 

{A,B}+ = AB + BA. 

Definition A. 2 A symplectic vector space (V, uS) is a finite dimensional real 
vector space V equipped with a distinguished bilinear form uj which is antisym- 
metric and nondegenerate, i.e. u>(u,v) = —oj(v,u) u, v G V and for every 
u 7^ G V there is a v E V satisfying u(u,v) 7^ 0. 

For example the standard space (JR 2N ,u>o) is a symplectic vector space with 
ujq{u, v) = v T au for all u, v G H 2N where a = ( 10)' 

Definition A. 3 A linear map S : V — > V on a symplectic vector space (V,u) 
is called symplectic or canonical if 

S*uj = oj 

where S*uj is the pullback 2-form given by S*lu(u,v) = u(Su, Sv) 

In the standard space (R 2Ar , loq) a matrix S is symplectic iff loq(u, v) = v T au = 
loq(Su,Sv) = (Sv) T aSu for all it, v G H 2N , that is S T aS = a. 
It follows immediately that det 2 S = 1 . To fix the sign we observe that the 
exterior power = u>o A ujq A ... A ujq with N factors is a volume form since loq 
is nondegenerate. We know that there is only one alternating iV— form on V 
which gives the volume spanned by the vectors u, with i = 1, ...N it is acting on. 



73 



74 



APPENDIX A. DEFINITIONS AND PROOFS 



This is the determinant itself or the determinant multiplied with an arbitrary 
nonzero number. 

If S is a matrix in H 2N then 

Q(ui, U2N) = S*Q(u±, U2N) = £l(Sui, SU2N) 

= c det[Sui, Su2n] = c det S detfiti, t*2jv] 
= det S Q(ui, U2n) 

All symplectic matrices therefore have determinant plus one. The proof was 
taken from Ref. [4]. 

Lemma A.l The trace of an N —mode Weyl operator W$, £ G R 2Ar is given 
by tr[W^\ = {2tt) n 5 2N {£), where 5 2N (£) = 5{£i) • ... • 6(£ 2 n) denotes the 
2 N— dimensional delta distribution. 

Proof : With the properties of the coherent states (see Section 3.1) we can 
calculate the trace of every operator in this overcomplete basis. 

tr[A] = ^ (n\A\n) = / dR(a)d9r(a) (n\a)(a\A\n) 

n=0 ^ n=0 ^ 

i r 00 

= - dft(a)d9(a) J2(a\A\n)(n\a) 

11 ^ n=0 

= i J dK(a)d9(a) (a\A\a) 

With the tensor product structure of the Weyl operators W% = ®jLi , 
£ = (r] 1 , . . . ,r] N ) T G ]R 2N ,rf G R 2 , it is sufficient to calculate only the one- 
mode case. Using the identity for the one mode W v = D{—0) for [3 = Vl ^ 2 

and furthermore D(/3)D(a) = e"^^*) D(/3+a), we start to calculate the trace 
of the Weyl operators. 

tr[W- v ] = ^ /"dSft(a)d»(a) (a\D((3)\a) 
1 



dK(a)d^(a) (0\D{-a)D(P)D(a 
- [ d#(a)d»(a) e -*W) (0|L>(/?)|0) 



7T 

= I e -W-3(/3) 2 J d 9(a) e ««0»)9W y dSR(a) e-^O^W 
= /" d»(a) e^^) | dft(a) e -^»«(«) 

= 2vr «5(\/2r/i) 2vr <f(\/2?72) 

7T 

= 27T 5(77l) <5(772) 

Immediately we see that for the N— mode Weyl operators irfVP^] 
tr[®? =1 Wj] = nf=i tr[W vJ ] = (2tt) n 11% as stated. 
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Lemma A. 2 (Baker-Campbell-Hausdorff) 

For M x M matrices or operators A and B 

e * Be -^B + l A,B ]+ ^; B K lAM f- m ... (A.1) 
For C = [A, B] and [A, C] = = [B, C] : 

e A + B = e A -e B -e-%. (A.2) 
Lemma A. 3 For every two positive N x N matrices A and B : tr\AB\ > 



Proof : For every positive matrix A the square root is denned by \[A := 0^/~DO T 
were D is the diagonal matrix of the (positive) eigenvalues of A and the 
orthogonal matrix which diagonalizes A. Since the trace is invariant under 
cyclic permutations it is tr[AB] = tr[\fA^/~AB\ = tr[\fAB\fA] > □ 

Lemma A. 4 The spectrum of a matrix C = AB with A and B square matrices 
and one of them invertible is the same as the the spectrum of the matrix D = 
BA. 



Proof : The eigenvalues of the product AB are determined by the equations 

= {AB - \it)xi = {A- X i B^ 1 )Bx i = B~ l (BA - X^Bxi. 

Since B~ l is nonsingular the same eigenvalue equations hold for the product 
BA. □ 



Lemma A. 5 The symplectic eigenvalues of a positive matrix are invariant un- 
der symplectic transformations S £ Sp(2N, R). 

Proof : The symplectic eigenvalues of a positive matrix A are given by the (usual) 
eigenvalues of —iaA, determined by [—iaA — Xit]xi = 0, with all positive A« 
symplectic eigenvalues of A. Let now S denote a symplectic transformation 
so that A = SA'S T . It follows that A' has the same eigenvalues as A since 
l-iaSA'S T - Xil]xi = S~ T [-iaA - Xil]S T Xi = 0. □ 
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